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PREFACE 


The G400PA mathematical formulations and resulting computer code described 
herein were developed by United Technologies Research Center (UTRC) under 
Contract NASl-16058, "Coupled Rotor Fuselage Vibration Analysis". This con- 
tract was through the Structures Laboratory of the USRTL (AVRADCOM) and 
Sikorsky Aircraft Division of the United Technologies Corporation was prime 
contractor. The initial development of the G400 analysis was conducted at 
UTRC by Dr. Richard L. Bielawa under Corporate sponsored independent research 
and development. Extensive refinements to the analysis were made under spon- 
sorship of the Langley Research Center of the National Aeronautics and Space 
Administration and the U.S. Army Air Mobility Research and Development 
Laboratory, Langley Directorate as part of Contract NASl-10960. This effort 
constitutes one of the contributing technologies to the Simplified Coupled 
Rotor Fuselage Vibration Analysis (SIMVIB) also funded through Contract 
NASl-16058. 
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Aeroelastlc Analysis for Helicopter Rotors 
with Blade Appended Pendulum Vibration Absorbers - 
Mathematical Derivations and ProRram User's Manual 


by 

Richard L. Bielawa 
United Technologies Research Center 


SUMMARY 


Mathematical development is presented for the expanded capabilities of 
the United Technologies Research Center (UTRC) G400 Rotor Aeroelastlc Analysis. 
This expanded analysis, G400PA, simulates the dynamics of teetered rotors, 
blade pendulum vibration absorbers and the higher harmonic excitations re- 
sulting from prescribed vibratory hub motions and higher harmonic blade pitch 
control. Formulations are also presented for calculating the rotor impedance 
matrix appropriate to these higher harmonic blade excitations. This Impedance 
matrix and the associated vibratory hub loads are intended as the rotor dynamic 
characteristics elements for use in the Simplified Coupled Rotor/Fuselage 
Vibration Analysis (SIMVIB) . 

Sections are included presenting updates to the development of the original 
G400 theory, and material appropriate to the user of the G400PA computer program. 
This material includes (1) a general description of the structuring of the 
G400PA FORTRAN coding, (2) a detailed description of the required input data 
and other useful information for successfully running the program, and (3) a 
detailed description of the output results. 


"k 

The research effort which led to the results in this report was financially 
supported by the Structures Laboratory, USARTL (AVRADCOM) . 





INTRODUCTION 


The G400PA helicopter rotor aeroelastic analysis described herein 
represents a generalization and extension of the original United Technologies 
Research Center (UTRC) G400 Rotor Aeroelastic Analysis which was first 
described in Reference 1. This analysis was originally developed in response 
to the unique aeroelastic characteristics of the composite bearingless rotor 
and represented an advancement in the state-of-the-art with regard to the 
modeling of rotors with time-variable, nonlinear structural twist and multiple 
structural redundancy. Since the publication of Reference 1, the G400 analysis 
has evolved into a family of analyses with a completely general range of 
applicability in rotor type (articulated, hingeless, teetered) and vehicle 
application (helicopters, propellers, wind turbines). 

The principal distinguishing features of the G400PA version of tlie analysis 

are: 

1. A capability for analyzing underslung teetered rotors. 

2. The implementation of the nonlinear equations of motion for two 
independent blade appended pendulum vibration absorbers. 

3. Inclusion of the dynamic and aerodynamic loads due to specified 
harmonic accelerations in arbitrary combinations of all six 
components of hub motion. 

4. The expansion of the blade pitch control description to include 
cosine and sine components of the (b-1) , b, and (b+1) harmonics 
of cyclic pitch angle, (where b is the number of blades). 

5. Provision for calculating impedance matrices, wherein each of the 
components of hub acceleration are systematically perturbed to 
obtain numerical partial derivatives of hub loads with respect 

to hub motion. 

Since the original G400 analysis is well documented in Reference 1, this 
report presents only the technical background, and descriptions and/or formu- 
lations pertinent to the G400PA features listed above. Additional sections 
are presented to facilitate the user of the actual (G400PA) computer program. 

One section gives a general description of the structuring of the computer 
program with descriptions of all the program elements. Other sections give 
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complete descriptions of all required inputs for the program as well as 
descriptions of the output results. Although these input and output 
descriptions overlap to a large degree those given in Reference 1, tliey are 
presented for completeness. The appendix updates, where appropriate, the 
material of Reference 1. 
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LIST OF SYMBOLS 


[A] 


Inertia coupling matrix for a single blade. 


[AA] 


The combined inertia coupling matrix for the degrees-of-f reedom 
of the two blades of a teetered rotor. 



n'th Fourier cosine and sine components, respectively, of inputted 
cyclic blade pitch control (higher harmonic control angles), deg. 




^66 


Submatrices comprising [A] . 


b 


Number of blades. 


c\j; , stj) 


Abbreviated notation for cosine and sine of respectively. 


C 

P 


Damper rate of explicit restraint of pendulum absorber about 
hinge, in-lb-sec/rad. 


DS1^,DS2^ ,DS3j^ , Integration constants used to calculate the blade internal loads 

DS4i^>DS5j immediately outboard of the flexbeam torque tube juncture, (ND) . 


e 


X2 coordinate of coincident flat-lag hinge or hingeless blade 
offset point, in. 


Function defining fraction of pendulum absorber mass to be included 
in mass of i'th blade segment, (ND) . 





Components of force in respective axis directions in the 5-coordina 
system, lb. 


h Distance rotor/blade apex point is below teeter hinge, in. 

us£ 




Identity matrix of dimension m. 


I 

P 


Mass moment of inertia of pendulum absorber about its flapping 
hinge, in-lb-sec . 


’"yio’^zio 


Mass radii of gyration of blade section about axes through and 
perpendicular to the spanwise (x^) axis and in the chordwise 
and thicknesswise directions, respectively, in.. 



LIST OF SYMBOLS (Cont'd) 


K Spring rate of explicit restraint of pendulum absorber about hinge, 

^ ft-lb/rad. 


K , K Effective spring rates defining the nonlinear elastic character- 

^1 ^2 istics of teeter stops, ft-lb/rad. 


I 


P 


m 


Distance from hinge to mass center of a pendulum absorber, in. 

2 2 

Blade mass distribution, lb-sec /ft . 


m 


•^A 


Reference blade mass distribution, defined by mass of the 3th 

2 2 

blade segment, Ib-scc /ft . 

Maas of pendulum absorber, Ib-sec^/ft. 


Components of blade moment about axes in the 5- coordinate system, 
arising from inertia, aerodynamic and/or concentrated mechanical 
effects, f t-lb. 

M^ The saturation moment used in the nonlinear elastic characteristics 

of teeter stops, ft-lb.. 

nj,nQ Blade segment numbers of inner and outer segments bracketing 

pendulum absorber mass. 

NEM,NFM,NTM Numbers of assumed flatwise, edgewise, and torsion natural 

"uncoupled" primitive modes, respectively. 


M 


X5-»y5.«.5 


Differential operator denoting dif ferentation with respect to 
time of quantities measured in the inertial frame, sec~^. 


Px5»Py5»Px5 


q 


Section shear load distributions in directions of axes in the 
5- coordinate system, Ib/in. 

General expression for a response variable deflection, (ND) . 


{q} 

‘^i 


Vector of blade degrees-of-freedom excluding flapping motion. 
Blade k'th edgewise modal response variable, (ND) . 

Blade i'th flatwise modal response variable, (ND) . 


'Ix5»'ly5,qz5 


Section moment load distributions about axes in the 5- coordinate 
system, In-lb/In. 
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LIST OF SYMBOLS (Cont'd) 




Ar, 


R 


S 

P 



^X3’^5*^Z5 


[ToUT^],[T2], 

[Tj] 

Av, Aw 


Blade j'th torsion modal response variable, (ND) . 

Blade spanwise coordinate, measured from offset, e, in 
direction, ft. 

i'th blade spanwise segment length, ft.' 

Rotor radius, ft. 

Amplitude of friction moment restraint of pencluluin absorber about 
hinge, ft-lb. 

Components of concentrated shear in directions of axes in the 
5- coordinate system, lb. 

Time, sec.. 

Coordinate system transformation matrixes relating rotating 
coordinate system deflection to the inertial frame. 

Deflection correction terms due to first order twist effects, (ND) . 


V^.Vy.V, 


Components of instantaneous lineal velocity of the nonrotating 
hub in the hub fixed coordinate system,, ft/sec. 


{V3} 


Vector of components of velocity of a point in the coned, lead- 
lagged, rotating coordinate system, ft/sec. 


AV, AW 


Deflection correction terms due to second order twist effects, (ND) . 


^5»y5>^5 


Components of the 5- coordinate system, defined to be rotating 
with the hub, but at the blade coned and lagged position, ft. 


{X3} 


Vector of components of displacement of a point in the coned , 
lead-lagged, rotating coordinate system. 


^ 10’ ^10 


Chordwise and thicknesswise position coordinate, respectively, 
of an arbitrary point within a blade section, in. 


’'“CG 


Chordwise distance of blade section mass center from the reference 
axis, in. 
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LIST OF SYMBOLS (Cont'd) 


Out of plane deflection in the unconed (and unflnpptjd) coordinate 
system, ft. 

Higher harmonic impedance matrix (partial derivatives of hub 
loads with respect to hub acceleration and HHC angles) . 

Total cone angle, consisting of built-in and flapping (or teetering) 
contributions, rad. 

Built-in blade precone, deg, rad. (as appropriate).. 

Flapping angle of blade appended pendulum vibration absorber, rad. 
Teeter angle (= 3-Bg) » rad. 

Angles defining the nonlinear elastic characteristics of teeter 
stops, deg. 

Deflection mode shape for the k’th edgewise normal mode, (ND) . 
Deflection mode shape for the I'th flatwise normal mode, (NU). 
Deflection mode shape for the j*th torsion normal mode, (ND) . 

Blade lead, angle, rad. 

Elastic torsion deflection angle, rad. 

Abbreviated torsion flexibility coefficients containing deflection 
dependency, rad/lb, rad/ft-lb. 

Blade pitch angle due to input control angle, rad. 

Built-in inclination of pendulum absorber hinge axis relative 
to local chord line, deg, rad (as appropriate). 

Total local blade pitch angle, rad. 

Total torsion deflection of flexbeam at junction point beyond 
built-in value, rad. 

Vector of excitations for twor-bladed teetering rotor degrees- 
of-freedom. 



LIST OF SYMBOLS (Cont'd) 



Component vectors of the E excitation vector relating to modal 
excitations and flapping excitation, respectively. 





X y 


,U) 


z 



The moment imparted to the teeter moment by the teeter stops. 
Blade azimuthal (angular) position, rad. 

(Nondimensional) uncoupled natural frequencies of i'th flatwise 
bending mode, k'th edgewise bending mode and j'th torsion mode, 
respectively, rad/sec. 

Components of instantaneous rotational velocity of the non- 
rotating hub in the hub fixed coordinate system, rad/sec. 

Rotor rotational frequency, rad/sec. 


SUBSCRIPTS . 


( ), 


( )t 


( ). 


( ) 


HHC 


( ) 


( )p.( ) 


PA 


( ) 


PAH 


( ) 


TS 


Structurally built-in parameter, or conditions of blade immediately 
outboard of juncture. 

Effect of dynamic origin. 

Due to elastic deformation. 

With respect to the nonrotating hub. 

Higher harmonic blade pitch control 

With respect to inertial space. 

Conditions at f lexbeam-torque tube juncture. 

Pendulum absorber. 

Pendulum absorber hinge. 

Pertaining to teeter stops. 
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LIST OF SYMBOLS (Cont'd) 


SUPERSCRIPTS 


Pertaining to the m'th blade of a teetered rotor 
Nondimensionalization by combinations of ni^^jR and/or 
Differentiation with respect to (Dt) . 

Differentiation with respect to (r/R) . 



SIMULATION OF UNDERSLUNG TEETERED ROTORS 


Basic Modeling Considerations 

For articulated and hlngelcHH rotor HystemH tlie Indlvldiia] rotor bladen 
are meclmnicaliy coupled to eucli other only througli motions of the huh. 
Consequently, these rotor types can be conveniently analyzed us single- 
bladed rotors whenever the assumption of Infinite hub rigidity is Involved. 

That is, without hub motion, the blades do not interact mechanically with 
each other. This is not true for teetered (two-bladed) rotors, however. 

Even for infinite impedance of the supporting shaft, each of the blades of 
a teetered rotor can. impart bending loads to the other through the root 
load restraint. Accordingly, teetered rotors must be analyzed, as a multiple 
(two) bladed rotor system. 

The method of implementing the teetered rotor capability in the G400PA 
analysis was selected based upon the already inherent use of the Galerkin 
method of solution. The selected method actually amounts to a "modified 
Galerkin" approach in that the assumed uncoupled bending modes used to 
define the blade elastic deformations actually satisfy only the geometric 
boundary conditions. More specifically, the modeling is based upon the 
following considerations: 

1. The rotor system is assumed to consist of a rigid body finite mass 
hub which teeters as a rigid body* and two blades which are each c antilever 
attached to this hub mass. Thus, the out-of-plane displacements for either 
blade consists of a rigid flapping part plus an elastic part based on 
"cantilevered" uncoupled mode shapes. For untwisted blades, this displace- 
ment distribution (for the n'th blade) assumes the following form (using 
the nomenclature of Reference 1): 

■' jr| *'i 

where y (r) are the uncoupled cantilevered flatwise bending mode shapes. 

In this Formulation, each of the modal response quantities are inde- 

pendent degrees-of-freedom. The advantage of this formulation is that symmetric 
bending rotor modes, which respond typically at even harmonics, are simulated 
exactly, and the antisymmetric rotor modes, which respond typically at odd 
harmonics, can be well approximated geometrically as "pinned" modes witli 
Equation (1). 
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2. The blade inplane deflections are simulated with no lead lag rotation, 
but with combinations of edgewise and flatwise bending as in a typical hinger 
less configuration. 

3. The mass of the hub, to which the blades are cantilever attached, is 
set to zero as a limiting process. The required moment balance about the teeter 
hinge (in the flapping direction) is achieved by properly combining the flapping 
equations for the two blades together with the moment imparted to the teeter 
motion excitation by the teeter stops. 

4. This formulation is inherently incapable of yielding correct bending 
moment calculations at the root region using mode deflection methods. This is 
because the formulation uses modes which do not satisfy the "load" boundary 
conditions at the root. Both the modified Galerkin formulation and the force-: 
integration method, however, are relatively insensitive to these boundary con- 
ditions. Consequently, the bending moments and stresses obtained using 
force-integration, should be reasonably accurate. 

The four considerations discussed above form the basis for the actual 
mathematical formulation discussed in the following subsections. In the first 
subsection, the load distributions resulting from the flapping motion of an 
underslung teetered blade are described. In the next subsection, the imple- 
mentation of the moment equilibration across the teeter hinge is presented. 
Finally, in the last subsection, the details of the implemerttation of the 
teeter stop load characteristics are presented. 


Aerodynamic and Dynamic Load Distributions 

The analysis presented in Reference 1, and implemented in the parent G400 
code, would be sufficient for formulating the teeter flap angle equation, but 
for the presence of the undersling distance, This distance is defined as 

the distance of the rotor-hub apex below the teeter axis. The approach taken 
for defining the additional loads due to undersling is to define an incremental 
position vector and to then form the appropriate incremental components of 
inertial velocity and acceleration, as measured in the blade coordinate system. 
Using Equation (35) of Reference 1 as a starting point and referring to Figure 1, 
one can write the Incremental (nondimensional) position vector as: 



( 2 ) 
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Figure 1. Schematic of Teetered Rotor System — Kinematic Definitions and 
Flapping Moment Equilibration 
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Using coordinate system transformations as described In Heference 1 and 
appropriate differentiation yields Llie following i.-xpresslon for Incremenla I 
velocity vector: 


« 

P 

0 

This vector then modifies the definition for blade section inflow angle 
used for defining the quasi-static angle-of-attack, and finally the blade air- 
load distribution. Further use of the coordinate system transformations, 
additional differentiation-and D'Alembert's principle, yields the following 
vector for the incremental (inertial) dynamic load distribution: 



^usl 





( I 

, 

1 A-(^-/3b) I 

K\ 

> = ~mh ^ 
. usC 



(4) 


Before the teetered blade root constraint can be implemented, the basic 
flapping equation must also be modified to account for the undersling distance. 
With reference to Figure 1, equilibration of the moment about the flapping hinge 
results in the following modified flapping equation: 


/' [(r-huse 


(5) 


'e 5 '5- 

u 

This modified equation is to be compared with Equation (31) of Reference 1. 
Note that p^^^ and p^^ are, respectively, the total load distributions in the 
x^- and z^- directions and contain the incremental loads given by Equation (4). 
The loads due to teeter stops are accounted for by the. general moment about 
the flapping hinge, this moment defined in greater detail in the following 

subsection. 
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Implementation of Teetered Blade Root Constraint 

As implemented by the G400 analysis for a single blade (infinite impedance 
hub), the blade equations of motion take the following form: 
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I A* 1 

In tliis form, J|he elasLlc mode acceleration subvector, iq f, and the flappln)’ 
acceleration, 3 , are linearly coupled, and all remaining linear and nonlinear 
terms in the various equations are grouped on the right-hand side as an excita- 
tion vector. . This vector is composed of a subvector for the elastic modes, 

1“ J, and the term for the flapping equations, This form of the equations 

is^appropriate for time-history solutions, but requires a simultaneous equations 
solution at each time step to decouple the accelerations. 

Before two blades are mechanically constrained by the teetered blade 
boundary conditions. Equation (6) must be generalized to the following 
augmented form: 


__qq 






q^ ' ^/3/3 



(7a) 


or, in a more abbreviated form as:. 


AA)rQ*} = {H} 


(7b) 


The teetered constaint can be conveniently implemented by first defining 
the following constraint relationship: 



The teeter constraint is applied to the equations of motion. Equation (7), 
with the following matrix operations: 
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(9) 


[t1[a a][t1{q*}= \t^\\3\ +|ah|^ 


The partitioning of Equation (9) is as follows: 


' A<'U I n 

^^q_i 2 - 

0 l-ASfi1 a|,2» 




( 10 ) 


Equation (10) is then in the same form as Equation (6) and the same 
procedure must be followed to achieve a time-history solution (decoupling of 
the accelerations and then numerical Integration). Generally, the simultaneous 
equations solution of Equation (10) is accomplished using partitioning methods; 
the details of this essentially mathematical operation are outside the scope 
of this report and are therefore omitted. 


Teeter Stop Moment Characteristics 

The moment applied to the teeter angle equation due to the teeter stops, 
^~TS’ essentially nonlinear function of teeter angle. This moment is 

characterized in the G400PA analysis by the five parameters defined graphically 
in Figure 2. For teeter angles between ± the teeter stop moment is, of 
course. Zero. When the teeter angle reaches the first limit stop angle, (i^, 
the limit stop imparts an initially elastic reaction with a spring rate of 
Kg^. Between this first limit stop angle and a saturation angle, 
reaction moment is assumed to be a generally nonlinear function of teeter 
angle and increases to the value of with a rate of Kg 2 » the saturation 
angle, At teeter angles above the reaction moment again obeys an 

elastic relationship. Generally, the saturation spring rate. Kg , will exceed 
the initial rate, Kg^. However, this implementation provides sufficient 
latitude for modeling other nonlinear relationships. Mathematically, this 
implementation can be defined as follows: 
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BLADE APPENDED PENDULUM VIBRATION ABSORBERS 


Basic Assumptions 

The inclusion of blade appended pendulum absorbers in the aeroraechanical 
description of the rotor blade requires the definition of additional equations 
of motion, as well as modifications to- existing ones. The following list of 
principal assumptions forms the basis of the GAOOPA modeling of blade appended 
pendulum absorbers: 

1. Provision must be made for simulating a maximum of two (2) Imlepiuulent 
pendular absorbers. State-of-the-art Installations of such vibration absorbers 
have included up to two absorbers, each tuned to a different rotor harmonic. 

2. Each pendular absorber consists of a specified hinge location, a rigid 
body mass (with c.g. offset and rotary inertia about that hinge), and a re- 
straint about the hinge accruing from a rotary spring, a rotary viscous damper, 
and a (nonlinear) friction damper (see Figure 3). 

3. The hinge location and orientation for each absorber is subject to 
the following constraints (again refer to Figure 3): 

a. The hinge axis must be within a blade section plane, e.g., 
perpendicular to the spanwise axis. 

b. The hinge axis may be both translated and inclined relative to 

the local chord line. This translation and inclination define 
the distance, angle, ©op^» shown in Figure 3. 

c. The hinge axis may be located at any spanwise station. 

4. The mass center of each absorber is nominally chordwlse balanced at 
the blade elastic (reference) axis, (i.e., no chordwise mass offset, FlOp^ “ 

in the undeflected position. Note that for nonzero values of inclination angle 
(of the hinge axis relative to the chord line), deflections of the pendulum will 
produce time variable chordwise c.g., offsets. 

5. To enable the analysis to make rapid parameter variations on the 
pendular parameters, the blade bending and torsion modes are calculated ex- 
cluding the mass of the pendula. Therefore, the effect of incremental pendulum 
mass must be explicitly accounted for in the blade dynamic equations. This is 
accomplished by adding the respective pendulum masses as Incremental mass 
distributions extended over the two blade segments which straddle that pendulum's 
radial location. 
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Figure 3. Pictorial of Typical Pendular Absorber Installation on Blade 
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6. Lastly, the pendula deflection angles are not restricted to "small" 
amplitudes. Thus, the realistic effects of amplitude saturation can be 
analyzed . 


Dynamic Loads 

The derivation of the dynamic loads acting on the pendulum follows the 
basic method described in the previous chapter wherein an incremental displace- 
ment vector is defined, and then with appropriate differentiation and coordinate 
transformations an inertial acceleration vector is formed. The displacement 
vector for a point mass hinged to an offset, inclined axis is given by: 


( COSPp 

* l"sin ®PAH sin/3p 
\ cos ®pAH 



( 12 ) 


where the total pitch angle of the hinge axis, 0p.„» is the sum of the blade 
pitch axis, 0, and the built-in pitch angle relative to the blade, QopA* 
resulting components of inertia loading are then given by the following equations: 


[-Sin /3p]^p - cos /3p+ 2 cos /9p ^p sin ©pA„ 

+ ()3+2 0) sin/SpCosSpAh] 

[-cos sin ©p;^„'^-2sin /9p^p 
+ sin )3p(l + 28+^p+5^)sin©p^„ 

-(zcos /3p(^p0+)S^p) + sin)3p(2^+V)) cos ©p^h] 

+ Z|Op^^[" 2 (Wg - C^+Zp) cos © 

+ {i+^^+ 2/3^ +2S)sin©]} 




-Mpa'i * 
PA ^o^rj 




^PA^ 


nnoA. 


(13a) 


(13b) 
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N 


-MpA^ 


PA 0 *^1 


{t [cos /3p COS 

-(2cos)3p(4pfi+/3^p)+sm'^pe) sin©p^„ -sin /3p(/3^+ 0+2/36) cos© p^h] 
+ Z|Op^j^[-(^^ + 2/9©+/9^)cos© -*(5sin©jJ. 


(13c) 


where the distribution function, is defined for only two of the blade 

spanwise segments, the centers form a straddling pair: the one closest In- 

board, r(n^), and the one closest outboard, r(n^), of the pendulum mass centiT 
location, r : (see Figure 4): 

A x\ 


^(np)-rpA 

r(np)-r(nj) 

^PA 

r(no)-r(rij) 


i i = nj 


; J:nr 


(14) 


Equations (13) and (14) are then in a general form appropriate for inclu- 
sion in the blade dynamic equations. Note that the nominal inertia load 
distributions, given by Equations (41a, b, c) and (44a, b, c) of Reference 
1 must also be modified since, at the two straddling radial stations, the 
mass distributions must be augmented by 


Equations of Motion 

The usual blade dynamic equations (those for modal responses, flapping 
and lead-lag motions) are easily modified by including the additional inertia 
load distributions as defined in the above subsection. The pendulum equations 
are formed by equilibrating the inertia moments about the hinge to the spring 
and damper restraints: 


TpiTp = e[cos/3pFzg-sini8^FxJ -Kp/3p-Cp^p-Sp sgn(^p) (15) 

where the concentrated loads acting on a pendulum mass center, and 
are given by: 





(16a) 
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Definitions and Moment Equilibration 





(16b) 



-sin 



+ COS 0pAH 



and where the quantities P c* sre components of tlie usual blade 

inertia loadings for the straddling segments modified by the following mul- 
tiplication: 


PpA - Z 




i=ni rnjjArjm 


P(ri) 


(17) 


The quantities, (Ap^^, Ap ^ and Ap are formed from Equations (13a, b and c), 
respectively, but with th^ factor ffjL/mQAr^) set equal to unity, and all the 
§ dependent terms deleted. All such g terms are instead grouped into the 
s?ngle inertia-acceleration term on the lift-hand side of Equation (15). Note 
that Ip is the effective rotary inertia of the pendulum about the hinge. The 
above development is the same for either pendulum and subscripting to denote 
pendulum index was omitted for clarity. 
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BLADE LOADS DUE TO HIGHER HARMONIC EXCITATIONS 


As originally formulated, the G400 aeroelastic analysis simulated a 
wide spectrum of "higher harmonic" blade loads. The term "higher harmonic" 
is used herein to denote frequency content in the range of blade number (b) 

X rotor rotation frequency (D) , and above. Indeed, because of the multi- 
harmonic nature of the local dynamic pressure of helicopter blade section in 
forward flight, higher harmonic airloads ensue even from the usual collective 
and simple harmonic cyclic blade pitch angles. Another explicit source of 
higher harmonic airloads is the considerable harmonic variability of the 
inflow velocities induced in the rotor plane by the vortex structure of the 
rotor itself. Reference (2) presents an in-depth description of this theory 
and its interactive implementation with the G400 analyses. 

Two further sources of higher harmonic excitation are considered herein. 
The first, prescribed harmonic motion of the hub , is required to evaluate the 
characteristics of the rotor coupling with the airframe. This coupling is 
accomplished in the Simplified Coupled Rotor/Airframe Vibration Analysis 
(SIMVIB) using the rotor impedance matrices computed in the G400PA code. 

The details of this calculation are described in a subsection below. The 
second source of higher harmonic excitation considered, higher harmon ic p itch 
control , is included to provide an important analytic capability for appli- 
cation to a timely area of helicopter vibration of research and development. 
Technical developments of each of these sources of higher harmonic excitation 
are presented in the subsections which follow. 

Prescribed Harmonic Hub Motions 

The required impedance matrix is the collection of partial derivatives of 
the six components of hub vibratory loads (3 shears and 3 moments) with respect 
to each of the six components of hub vibratory acceleration (3 translations and 
3 rotations). This calculation is generally implemented by first achieving 
harmonic responses with selective perturbations in each of the components of 
acceleration and then forming numerical partial derivatives of the hub loads. 
The actual hub loads are calculated using the force-integration method (see 
Reference 3) together with the total resulting inertia, gravity and aerodynamic 
blade load distributions. The details of this calculation do not represent 
new technology and, consequently, are omitted herein. 



Of more importance are the calculations of the incremental airloads and 
dynamic loads due to harmonic hub motion. For this purpose, the components 
of inertia velocity and acceleration, in the blade coordinate system are 
required. Generally, the hub degrees-of-freedom are defined in the non- 
rotating (with rotor speed, fi) body-fixed axis system; therefore, the basic 
degrees-of-freedom are the instantaneous velocity components along (or about) 
coordinate axes fixed to the "nonrotating", but generally moving hub. Because 
of this definition of basic hub degrees-of-freedom, the actual deflections 
(translations and rotations) are subjtjct to the order of Integration and tlie 
appropriate Euler angles must be suitably defined and Integrated in this 
context. Furthermore, recourse must be made to tlie Theorem of Coriolis 
(Reference 4) to obtain the correct differentiation with respect to time 
as measured in an inertial space. The material to follow outlines the 
approach for making these differentiations and obtaining their components 
in the blade coordinate system. Much intermediate calculation is omitted 
for clarity, and the final results include the previously considered effects 
of undersling (teetered rotor feature) and thickness noise mass center off- 
sets (pendulum absorber feature) . 

Consider the vector sum for a point mass at a blade section: 

Xj--Xo + X, + X2 (18). 

where is the position vector with respect to inertial space. X^ is the 
position vector of the nonrotating (with respect to fi) hub: 

(19, 

[T ] is' the coordinate transformation due to hub orientation motion and^X, I 
o * h I 

is the position vector of the hub in the hub fixed system. Consistent with 

the definition of the hub centered coordinate system, [T ] and j X, 1 can be 

o * h • 

considered to be unity and null, respectively, at any instant in time, but 
have finite time derivatives, position vector of the offset/ 

blade root point: 

W=[T„][T,]{e} (20) 

where [T ] is the coordinate system transformation due to rotor rotation, and 
|e I is the position vector from the hub center to the blade root point. 1^2^ 
is the position vector due to blade element displacement within the blade 
coordinate system: 

( 21 , 
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whero [T^] is tht: coordinate transformation due to fiai)pinj> and iead-iaj; alxnil 
a common point, and tlie position vector of the point mass within tlu‘ 

blade "5" coordinate system. Specific definitions for the coordinate trans- 
formations and position vectors are as follows: 


Po] = 


0 “‘*'2 

0 “w 


— CUvi 


Wy 


(22a) 


w= 


C\|/ -stj/ 0 

Srp 0 

0 0 I 


(22b) 




{e}: 


I -S -/3 
S I 0 
/3 0 I 


0_ 

“husJ 


(22c) 


(22d) 


{^ 5 }=' y,oC® -Z,oS® +€y.l 

yioS® -"ZioC® 


(22e) 


where 6 , € and € are the elastic defledtions in the x^, and z. directions, 
X y z 3 j j 

respectively, and are functions of the modal responses, structured twist, etc. 

Furthermore, for convenience, a matrix [T ] is defined which transforms from 

the nonrotating blade root to the blade "5" coordinate system: 




(23) 


Other operations on the matrix which are useful: 


r ,-i r I 8 ^ircv// 0 l 

M 5 I = -8 1 0 -S</^Cv#^0 

0 'Jl 0 'J 


(24) 
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0 
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( 25 ) 


The total inertial velocity vector, is defined using the Theorem of 
Coriolis and is given in matrix form as follows: 

^h][T|]{^}- ^[[^o]h] * k][y {X5} 


(26) 


where the operator p denotes differentiation with respect to time, as measured 
in the inertial coordinate frame. 

The terms underlined are deleted because they do not contribute to hub 
motion. Note that in this formulation, theactual hub deflections (transla- 
tions and rotations) are taken as zero at any instant in time. 

Transformation into components in the "5" coordinate system is accomplished 
by premultiplication by [T^] , and retention of only those terms relating to 

hub motion and/or undersling results in the following: 




Vx ( c 8 S ) + Vy ( S :// + 8 C <j/ ) +/3 V2 

+ Sc»(/)+Vy(c»//- 8s>|/) 
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Equation (27) then gives the components of incremental velocity at the local 

blade section with which the tangential and normal components or air velocity, 

and U , respectively, can be |ppropriately modified. The total blade 

section £n flow angle, 4 > (= tan (U /U )), is then formed in the usual manner 

p T 

and included in the section angle-of-attack. In this manner, the effects of 
prescribed hub motion are included in the airload distribution. 

The inertia (dynamic) load distributions require the second time deriva- 
tive with respect to inertial space. The Theorem of Coriolis is now applied 
to the velocity vector, Equation (26), rather than the displacement vector: 




+ 2 
■ + 
+ 2 


[To][T|] 4To][f,] {S}^ [To][Tj{y,} 


. (28) 


The underlined terms are deleted because they do not contain contributions 
accruing from hub motion. After the indicated operations have been performed, 
the following expressions for incremental components of inertia load distri- 
bution are obtained: 
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In both Equations (27) and (28) the expressions for hub velocity and 
acceleration are completely general. For present purposes, these components 
of velocity and acceleration are to be prescribed as simi)le (but lilglier) 
harmonic functions: 


I S(n»^) (30) 

where j = x, y, or z, and where n is equal to a prescribed multiple of the 
number of blades.. To minimize the computational requirements for calculating 
these ^rigonometric functions, the following identities are utilized: 

1 _£!_ c(m</) 

1 Cff/ 

2 l-2(S«//)^ 

3 4(Ct//)^-3(Ci//) 
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Higher Harmonic Blade Pitch ControJ 
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The implementation of higher harmonic pitch control (HHC) in the G400PA 
analysis is a relatively straightforward task. For this type of higher 
harmonic excitation, the incremental blade pitch angle is assumed to be of 
the following form: 

■An, C(bi^/)-Bn5S(b.//) (32) 

-An+i.C(b+l).i/-Bn+|S(b+l)'<' 

s « 

The amplitudes in this partial Fourier series are generally prescribed, 
or are perturbed in order to obtain the partial derivative matrix of hub load 
with respect to HHC pitch amplitudes. The incremental pitch control angle 
given by Equation (32) is then ^scd to^^efine the blade airloads and to com- 
lilete the calculations for 0 , 0 and 0 used In forming the blade Inertia 
load distributions. 



Impedance Matrix Calculation 

The impedance matrix is defined as the array of partial derivatives of 
components of hub loads with respect to components of higher harmonic excita- 
tion variables. 

I<i<6 (11) 

l<j<9 

These X- and Xc terms are, respectively, the cosine and sine components of each 
'"i ®i 

of the SIX components of hub load, X^, and similarly qc . and q^ . are, respect- 
ively, the cosine and sine components of the higher harmonic excitation vari- 
ables, q.. Specifically, these consist of each of the six components of hub 
vibrator^ amplitude (cosine and sine) and three higher harmonic control pitch' 
amplitudes (cosine and sine). The calculations for these derivatives are 
accomplished numerically by first obtaining responses which are converged 
to periodicity without, and then, with perturbations in each of the higher 
harmonic excitation variables, q.: 


f zl- 

d(Xc. ,X,.) 


d(qc..Q-.) 

L *j J 


aXj _ Xi(qj ■fAqi)-Xj (q)) 
^ Aq, 


(3A) 


In practice, the perturbation amplitudes are selected to achieve calcula- 
tions of the partial derivatives which approach limit values. Due to round-off 
characteristics of the computer and the inability to achieve exactly converged 
responses, there is generally a practical limit to how small the various Aq . 
values can be. Some care and experimentation are therefore warranted in th^ 
selection of the Aq. perturbation levels. 
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GENERAL DESCRIPTION OF THE G400PA PROGRAM STRUCTURE 


The G400PA Rotor Aeroelastics Analysis is an extended, special purpose 
version of the parent G400 analysis described in Reference 1. The G400 analysis 
was revised to include the dynamics of blade appended pendular vibration 
absorbers, and the additional blade loadings (both aerodynamic and dynamic) 
due to prescribed sinusoidal hub motion. This version includes impedance 
matrix calculations, teetered rotor simulation, spanwise variable airfoil 
data, direct coupling with the UTRC Rotor Inflow Anlaysis (UTRC Program F389) 

and the use of higher harmonic pitch control angles. The eigensolution and 

transient spectral stability analysis features are deleted in this version 
of G400. The program was coded in FORTRAN IV and developed on the UNIVAC 
1110/81A computer system and has been converted to the 1BM370 and Langley 
Research Center CDC systems. 

Program Structure 

The G400PA program is structured in a generally conventional hierarchical 
fashion with an initial routine (MAIN), appropriate modular levels and a col- 
lection of utility subroutines and/or functions. The subroutine calling 
structure is presented in the tabular outline given below; this functional 
structuring also serves to list the elements (subroutines and/or functions), 

contained in G400PA. Note that the hierarchial structure is alternately indic- 

ated using indentations with indicating symbols and/or nested parentheses: 

Note also that this table lists the subroutines in alphabetical order within 
any given hierarchical level and not with respect to calling order, size or 
importance. The subsection which follows Table I presents a brief description 
of each of the G400PA elements. 
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TABLE I 


HIERARCHICAL LISTING OF G400PA ELEMENTS 


• MAIN 
■ NIAM 

ABNDGTT (GJR, QUADFT) 

A INISH2 

I ADMCOF 
0 DPCHEK 

0 PABSSU (ADMCOF) 

A INTCOF (FACT12, LNAERO) 
A LNAERO (BLIN5, SHLDM) 

A LOAD 
A MODEIN 
A MYKSTD (GJR) 

A PRNT 
A PTFLLG 

A RESETQ (QQPSET) 

A SPNWIZ 
A SPRINT 
A TRFLXl 
A TWOF 
■ MOTION 
A MAJITR 

0 F389SU 

▼ BLIN5 

▼ INTERP 

▼ LNAERO 

▼ SHLDM 

▼ STALN 

0 QSTHRM (HARM) 

0 RESETQ (QQPSET) 
0 TMARCH 
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TABLE I (Cont'd) 


T ADAMS 

▼ AERPRF (ALFCOR (HEYSDL(HEYSK) ) ) 
T BLADEL 

FACT12 
^ HBINRT 
HHCANG 
^ QQPSET 
.^RDHANL 
SETUP 
^ SPANS 

I BLIN5 ■ 

I UMl’RTB 
I HYSDMP 
I SHLDM 
I UNCOEF 

f BLIN5 
• ^ LNAERO 
f SHLDM 
f STALM 
f STALN 
IT ULSTAL 

^ TABLE 
T GROUTS 

▼ EXTRAP 

▼ QPPTST 

▼ STRSS 

T TEETR (GJR) 
f TRMRLB (GJR) 


A HUB SUM 



Brief Description of Program Elements 

ADAMS 

This function implements an algorithm for time integration of any of the 
time dependent variables. , 

ADMCOF 

This subroutine evaluates the coefficients required for the time integra- 
tion algorithm. 

AERPRF 

This subroutine completes the calculations for and outputs the summary 
of aerodynamic performance quantities. 

ALFCOR 

This subroutine calculates the corrected rotor angle-of-attack due to the 
proximity of a ground plane or wind tunnel walls using the theory of Reference 5 

BLADEL 

This subroutine provides ■ the computational loop over- number of blades in 
forming the blade response equations. The upper limit on this hoop is (1,2) 
depending on whether the rotor is of (unteetered, teetered) type. 

BLIN5 


This subroutine does a tri-variant table loop up. of the airfoil section 
coefficients. The three variables are angle-of-attack, Mach number and 
spanwise station. 

BNDGTT 

This subroutine calculates the flatwise and edgewise bending stiffness 
characteristic of the torque tube component of a bearingless rotor, as 
required for the redundant analysis calculation (see RDANAL) . 

GROUTS 

This subroutine is a compact simultaneous equations solver used for non- 
teetered rotor configurations. It uses the Grout Reduction method described 
in Reference 6. 
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DMPRTB 


This subroutine does a table look-up to implement the nonlinear lag damper 
angular rate vs damper moment characteristics. 

DPCHEK 

This subroutine tests the input integration steps size for acceptable 
accuracy, and. automatically decreases it if the value is too large. 

EXTRAP 

This subroutine effects a "static" solution on any degree-of-f reedom 
whose natural frequency is sufficiently high to approximate the response 
neglecting the twice time differentiated term in that degree-of-f reedom* s 
governing equation. 

FACT12 

This subroutine performs an interpolation function calculation in the use 
of input tabulated variable inflow. 

F389SU 

. This subroutine sets up the various blocks of data needed by the (UTRC F38y) 
variable inflow branch of the coupled rotor/fuselage vibration analysis and 
writes them to appropriate data files. A detailed description of the approp- 
riate data transmitted between the G400PA and F389 portions of the analysis Is 
given in Reference 2. 

GJR 

This utility subroutine optionally obtains simultaneous equation solutions 
and/or matrix inversions using the Gauss-Jordan Reduction method. 

HARM 

This utility subroutine performs a Fourier (harmonic) Analysis of any 
time history string of data. This harmonic analysis uses a negative cosine 
and sine definition for the harmonic components. 

HBINRT 

This subroutine calculates components of inertial acceleration, as measured 
in the rotating blade coordinate system, from the prescribed harmonics of hub 
motion. 



HEYSDL 


This subroutine calculates Interference factors needed for subroutine 
ALFCOR. 

HEYSK 


This subroutine calculates induced velocity functions for subroutine llEYSDh 
HHCANG 

This subroutine calculates contributions to the total blade pitch angle 
and its first and second time derivatives accruing from higher harmonic control 
angle inputs. 

HUBSUM 

This subroutine provides a printout of the rotor hub loads and impedance 
matrix calculations which are transmitted to the coupled matrix/fuselage vibra- 
tion analysis program. 

HYSDMP 

This subroutine calculates the increment to blade edgewise bending moment 
to account for hysteretic structural damping. This formation of structural 
damping is dependent on edgewise deflection and the signs of rate and accelera- 
tion, but not their magnitudes. 

INISH2 

This subroutine initializes arrays and logic variables, and nondimension- 
alizes parameters, as required for the time-history solution. | 

INTCOF 

This subroutine forms the modal integration coefficients. 

INTER? 

This utility subroutine provides linear interpolation calculations for 
subroutine F389SU. 

LNAERO 

This subroutine calculates partial derivatives of airfoil section co- 
efficients with respect to angle-of-attack and Mach number. 
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MAIN 

This is the main program element. Its explicit functions are to input . 
the tabulated airfoil section data and to call subroutines NIAM and MOTION. 

MAJITR 

This subroutine provides the structuring needed to produce a trimmed 
configuration, and to calculate the impedance matrices. Specifically, 
this subroutine provides the iterative loop structure to perturb various 
system parameters in a systematic manner to form numerical partial deriya- 
tives, both those needed for trim and those needed for the impedance matrices. 

MODEIN 

This utility subroutine inputs the blade bending and torsion mode shapes 
and their derivatives with respect to span. 

MOTION 

This subroutine controls the time-history solution flow. 

MYKSTD 

This subroutine calculates the static flatwise bending deflection and 
spanwise derivative distributions for the blade for unit load on the inboard 
and of the cantilever mounted torque tube (for bearingless rotor rotor applica- 
tions only) . 

NIAM 

This subroutine performs the following functions: 

1. completes the input of the required data; specifically (a) the 
inertia, elastic, geometric and other operational data, (b) the blade 
mode shape data, (c) the variable inflow data, and (d) the saved 
initial conditions, from an appropriate data file. 

2. performs some initializations and/or nondimensionalizations of 
logic and system parameters. 

3. calculates the deflection correction arrays which accrue from 
structural twist. 
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PABSSU 


This subroutine sets up parameters and provides appropriate nondimensionu 1 
izations as required for the prescribed hub motions and simulations of the 
pendular absorbers. 

PRNT 


This subroutine provides an echo print output of the Part II input data 
which pertains to the Inertia, Elastic, Geometric and other Operational 
Data. 

PTFLLG 

This subroutine calculates the parameters used to define the pitch-flat, 
pitch-edge and pitch-flap couplings. 

QPPTST 

This subroutine tests the system degrees-of-freedom for numerical instab- 
ilities. The criteria used to identify such an instability are the occurrence 
of three sign changes of increasing amplitude in three time steps. 

QQPSET 

This subroutine sets the blade modal responses and modal rates for a 
given blade from the vectors of system response amplitudes and rates, 
respectively. 

QSTHPM 

This subroutine performs harmonic analyses (using subroutine HARM) of the 
blade modal responses, pendular absorber responses, hub shears and moments, 
and blade stresses, after the responses have converged to periodicity. 

QUADFT 

This subroutine fits a quadratic function to any set of three ordinates 
corresponding to three equally spaced abscissae. 

RDNANL 

This subroutine makes the redundant anlaysis calculations needed for 
simulations of bearingless rotor configurations. 
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RESETQ . 

This subroutine places selected terminal conditions into an array and 
writes them to a data file for use as initial conditions for subsequent runs. 

SETUP 

This subroutine sets up the inertia coupling matrix and excitation vector 
for the modal responses of a given blade. , 

SHLDM . 


This subroutine analytically approximates the NACA 0012 airfoil aero- 
dynamic coefficients as functions' of angle-of-attack and Mach number. 

SPANS 

This subroutine implements the loop over spanwise station in forming 
the blade response equations. The spanwise inertia, aerodynamic and gravity 
load distributions are formed in this subroutine. 

SPNWIZ . 

This function performs a numerical integration between blade section 
centers of a specific integrand type as required for forming the deflection 
corrections functions due to structural twist. 

SPRINT 

This subroutine outputs (as optionally requested) the spanwise integra- 
tion coefficients. Although most of these coefficients were required only 
for the (deleted) elgensolution, some are used in the time-history solution. 
For completeness all coefficients were retained. 

STALM 

This function approximates the variation of static moment coefficient 
stall angle with Mach number for the NACA 0012 airfoil. 

STALN 


This function approximates the 
to close approximation) coefficient 
NACA 0012 airfoil. 


variation of static normal force 
stall angle with Mach number for 


(lift, 

the 
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STRSSS 


This subroutine calculates the spanwise stresses and Intesrateil hub loads 
optionally using the force-integration or mode deflection metliods. 

TEETR 

This subroutine combines the inertia coupling matrices and excitation 
vectors for a two-bladed rotor to effect the teetered rotor boundary condi- 
tions at the blade roots. 

TMARCH 

This subroutine controls the solution flow for obtaining the time-history 
solutions. It furthermore tests for numerical instabilities and convergence 
to periodicity. 

TRFLXl . 

This subroutine calculates the torsion deflection characteristics of the 
flexbeams (for bearingless rotor configurations) due to spanwise varying 
torques resulting from the concentrated shears and moments at junction point 
of the flexbeam, torque tube and blade proper. 

TRMRLB 

This subroutine calculates the final partial derivative matrices used 
to obtain a trimmed flight condition. 

TWOF 

This function performs a least-square curve fit calculation on blade 
twist to facilitate subroutine NIAM in the numerical differentiation of 
blade twist. 

ULSTAL 

This subroutine performs a table lookup of the input airfoil data to 
obtain (1) derivatives of section unstalled lift and moment coefficients 
with respect to angle-of-attack, and (2) static life and moment stall 
angles, for use in subroutine UNCOEF. 

UNCOEF 

This subroutine implements the UTRC synthesized, unsteady stalled air- 
loads formulations, as described in Reference 7. 



PROGRAM INPUT DESCRIPTION 


The required input to the program consists of the following major data 
blocks in order of loading: 

I. Airfoil Data 

II. Inertia, Elastic, Geometric and other Operational Data (in Loader 
Format) 

III. Blade Mode Shape Data 

IV. Harmonics of Variable Inflow 

Details for preparing the data for each of these blocks are given in the 
subsections which follow. An additional subsection provides information for , 
facilitating program operation and improving its efficiency. 


1. Airfoil Data 

This data block generally consists of one or more sets of tables of two- 
dimensional lift, drag and pitching moment coefficients versus angle-of-attack 
for various Mach numbers. for up to four (4) arbitrary spanwise locations. 
Additionally, if • unsteady aerodynamics are used, the static stall angles 
and linear coefficient slopes for both lift and pitching moment are included 
in this table. 

Input Format for First Card(s) 

While actual set-up of this data block follows a basic format (described 
below), specific variations are required on the first card(s) of this block 
depending on optional usage. These variations denote whether multiple data 
sets are to be input for respective spanwise locations, a single set is to 
be input for use on all spanwise locations, or an analytic representation 
of the NACA 0012 airfoil is to be used for all spanwise locations. Each of 
these optional usages are described below. 

Mu]^t^l^e_S£aiiw^s£ Aijr f^ij.s 

For the case of distinct airfoil characteristics being defined at up to 
four (4) multiple spanwise locations, the first card image format is as follows 



card //lA: 


f 


NZ NRCL NRCD NRCM TlTLE(optional) 


(412, A72) 


where NZ would be ignored so long as the sum of the quantities NRCL, NRCD 
and NRCM is 4 or greater. The quantities NRCL, NRCD and NRCM are, respectively,, 
the number of radial stations for which multiple c^, c^ and '^mc /4 ^^rfoil data 
are to be inputted, as appropriate, each with a minimum value. of 1 and a maximum 
value of 4. Note that at least one of these three inputs must have a value of 
2 or, greater.. For multiple spanwise section properties, an additional card, 
following the one described above, is then required, which begins the input of 

the c„ airfoil data: 

t . . 


card //IB: 


where NZ(1) is the number of Mach number Cor which groups of c^ data are to 
be read in for the first radial station; RADCL(l) is the nondiitiensional radial 
station at which the airfoil data is defined. 

^ingl^e_A^r^ol^l_D^S£r^_ti£n_ 

For the case of a single airfoil to be used for all spanwise locations, 
a single first card image is input. This card is similar to the card //lA 
described above, except the quantities NRCL, NRCD and NRCM are input as zero 
(or the columns are left blank): 


card. //I; / NZ(1) 0 0 0 TITLE(optlonal) (412, A72) 


In this case the card is interpreted as the first card of the c data 
with the RADCL(l) information omitted (see description above for card //IB) . 

An£l^tJLc_Aj.rf^o^l_D£S£r^p_ti^n_ 

For those optional cases wherein the analytic NACA 0012 airfoil option 
is invoked (see location 996 discussed in the following subsection) the first 
card image must be a single card with blank or zeroed columns 1 through 8. 

For this option, the remainder of the airfoil data block of data is omitted. 


NZ(1) RADCL(l) TITLE(optional) 


(12, F8.0, A70) 
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Input Format for Subsequent Cards 

For those cases wherein tabulated airfoil data are to be input, the 
subsequent inputs continue the airfoil c data initiated with the itlA and/or 
//IB cards. Thus, the card image set-up defined below is then input NRCL 
times (or only once, if NRCL = 0): 


card //2: 

/ J N M A(l) CL(1) A(2) CL(2) ...A(4) CL(4) 

(12, 10F7.0) 




cards //2+: 

/ A(5) CL(5) ... 

(F9.0, 9F7.0) 


• 

• 

• 



/ . . .A(N) CL(N) ALSTAL DCLDAO 

- 

(F9.0, 9F7.0) 


where: J is the number of data entries to be input into each Macli number group 

N is the number of angle-of-attack -c abscissae-ordinate pairs to be input. 
Normally a maximum of thirty-four (34) a pairs may be input; thirty-three 
(33) pairs are input if the unsteady option is chosen. M is the Mach number 
appropriate to the data group. A(i) are the N angle-of-attack abscissae in 
degrees and CL(i) are the N lift coefficient ordinates. ALSTAL and DCLDAO 
are, respectively, the static stall angle, in degrees, and the lift curve 
slope at zero angle-of-attack, in per degree; these items are needed only 
if the unsteady airloads option is invoked. Note that J is fixed-point 
formatted, but N and M are floating point formatted. 

Cards 2. and 2+ are repeated for each successively higher Mach number. 

A maximum of 12 Mach numbers is allowed and the lowest and highest Mach 
numbers need not define the total working range as the search technique uses 
the boundary data for Mach numbers beyond the range input. Thus, repeated 
data for zero and supersonic Mach numbers are not needed. The lowest Mach 
number input must contain an angle-of-attack range from -180° to 180° or from 
0° to 180° depending on whether or not unsymmetric airfoil data is being input; 
all higher Mach number data need extend only from -30° to 30° or from 0° to 
30° in a similar manner. 





The general format described above is repeated for the and 
subblocks in that order but with either card image //IB or //lA used, depending 
on whether multiple airfoil section data are or are not being input (and 
used). The static stall angles and aerodynamic coefficient curve slopes at 
zero angle-of-attack are deleted for the c^j subblock. Lastly, note that the 
total storage allocated for the combined airfoil data is AOOO (which is 
somewhat less than 2000 abscissa - ordinate point pairs). The maximum storage 
should be adequate for most applications and can be "budgeted" among the 
various spanwise sections, as appropriate. 


II. Inertia, Elastic, Geometric 
and Other Operational Data 

This data block includes those items used to define the more detailed 
dynamic features and/or those which.' are most likely to vary from case to 
case. The card image format for these data is as follows: 


NN L DATA(L) DATA(L+1) ... DATA(L+A) (12, lA, 5F12.0) 


where: NN is the number of data items to be input on a given card; NN must 
not exceed 5. L is the location or identifying number of the first data item 
on the card columns 3-6 right adjusted. DATA(L+1) represents the various 
data items on the card image, columns 7-18, 19-30, 31-A2, A3-5A, and 55-66, 
in the floating point format. The locations or identifying numbers for the 
various data and control items are listed below along with definitions and 
other pertinent comments; note that some data locations are intentionally left 
blank. 



Locat ion Item 

1 'i2R 

2 R 

3 P 

A a 

CO 

5 b 


De scription 

Rotor tip speed, ft/sec. 
Rotor radius, ft. 

2 A 

Air density, lb-sec /ft . 
Speed of sound, ft/sec. 
Number of blades. 


A6 




Description 


Nondimensional offset distance of start of 
deformable and/or deflectable portion of 
rotor blade (hinge offset for articulated 
rotors), e/R. 

Tip loss, used to define equivalent momentum 
area pnd three-dimensionality corrections 
to computed two-dimensional airloads near 
the blade tip. 

Number of blade segments used to define span- 
wise variable arrays. 

Azimuth increment used in the numerical 
integration of the dynamic equations, deg. 

(See section on general information for 
efficient program usage.) 

Number of "flap trials", i.e., maximum number 
of rotor revolutions for which the blade time- 
history will be computed In an attempt to 
obtain convergence to periodicity. if a 
transient response is desired for only a 
portion of one rotor revolution the program 
will compute a time-history solution for any 
nonzero fractional N value input. An identic 

r 

ally zero value will cause the time-history 
solution to be by-passed entirely. 

Flapping tolerance to within which the aero- 
elastic/dynamic responses must repeat on 
successive revolutions in order for the 
motion to be considered converged to period- 
icity. The tolerance applied to lead-lag 
motion is equal to 5c^. 

Longitudinal cyclic pitch, coefficient of 
minus cosi^ term in Fourier expansion of 
blade control pitch angle, deg. 

Lateral cyclic pitch, coefficient of minus 
sin<j) term in Fourier expansion of blade 
control pitch angle, deg. 

Blade collective pitch angle as defined at 
the 75% radius, deg. 



Location 

15 

16 

17 

18 


19 

20 

21 

22 

23 


Item 


X 


V 

0 


(Control) 


N 

cut-out 


d cut-out 




print 


Definition 

Mean rotor inflow ratio. 

Forward flight velocity, kts. 

Rotor area solidity (bc/nR) . 

Increment added to all values of c^ obtained 
from tabulated airfoil data or from the 
analytic NACA 0012 data. Airfoil data 
generally correspond to smooth wind tunnel 
models and Ac^j is often used to adjust for 
the higher drag of production blades; a 
commonly used value of Ac^j^ is 0.002. 

Make greater than zero (1.) for first case 
or when new blade modal data are to.be input. 
Program automatically sets this control number 
to (-1.) after each loading of modal data. 

Number of blade segments, starting at inboard 
end and defining the cut-out region, for which 
the lift and moment coefficients are set to 
zero. 

The drag coefficient used on the first 

N segments, 

cut-out 

Built-in linear blade twist angle; i.e., dif- 
ference between tip and root built-in angles, 
positive when tip angle is greater (L.D. up) 
than root angle, deg. 

Azimuth increment used to present printed 
output of various pertinent aerodynamic, 
dynamic and elastic load distributions as 
well as aeroelastic responses and stresses, 
deg. This input quantity should be an in- 
tegral multiple of location 9, deg. 



Location 

Item 

Description 

24 

% 

"Direct" value of pitch-flap coupling (AQ/A8) . 
A nonzero value will suppress a calculation 
of this quantity from the input pitch-horn/ 
push-rod geometry. 

25 

^6 

Pitch-lag coupling (AO/AiS). 

26 • 

^LD 

Viscous lag damper coefficient, ft-lb-sec. 

27 

c 

Blade chord if chord is constant, otherwise 
omit, ft: 

28 

NFM 

Number of flatwise bending' modes to be used 
(4 max) . 

29 

NEM 

. Number of edgewise bending modes to be used 
(3 max). 

30 

NTM 

Number of elastic torsion modes to be used 
(3 max) note that the total number of dynamic 
degrees-of-freedom is limited to 10. Thus, 
if the articulation degrees-of-freedom, 8 
• and/or 6, are used, or if the rigid feathering 
d.o.f. is invoked, NFM, NEM and NTM will be 
automatically changed to keep the totaJ degree 
of-freedom to no more than 10.). 

31-34 

“wi 

Flatwise modal frequencies, nondimensional 
with respect to 0, in ascending modal order. 

37-39 

(0 

Edgewise modal frequencies, nondimensional! 

with respect to fi, in ascending modal order. 

■ 

40-42 

<^0. 

J 

Torsion modal frequencies, nondimensional 
with respect to fl, in ascending modal order. 

43 


Second harmonic cyclic pitch coefficient of 
minus cos 2ijf in Fourier expansion of blade 
control pitch, deg. 

44 

1 


Second harmonic cyclic pitch coefficient of 
minus sin 2tj) in Fourier expansion of blade 


control pitch, deg. 



Item 


Description 


Location 

45 g 

46 (Control) 

47 (Control) . 

48 (Control) 


49 


(GJ) 


root 


50 


K 


Qroot 


51 (Control) 


Acceleration due to gravity, ft/sec ; a nega- 
tive value implies inverted flight. 

Make nonzero (1.) if airfoil data for a 
nonsymmetric airfoil are to be used. 

Make nonzero (1.) if the total ( transient ) 
time -history is to be output; i.e., responses 
calculated before convergence to periodicity 
is obtained. 

Make nonzero (1.) if the modal responses and 
hub shears an d moments are to be (negative) 
Fourier analyzed after periodicity has been 
obtained. 

Torsional rigidity at the blade root, Ib/ft^. 

Torsional spring rate connecting root of blade 
to fixed structure to represent control system 
flexibility, ft-lb/rad. A nonzero value will 
automatically introduce the rigid-body feather- 
ing degree-of-freedom as an addition "torsion 
mode". Note that this "mode" will inherently 
couple with the NTM normal elastic torsion modes 
at frequencies both below and above the input 
torsion frequencies (locations 40-42); hence, 
a smaller integration interval will generally 
be required. 

Make nonzero (1.) to simplify the numerical 
spanwise integration techniques from the 
nominal trapezoidal rule to rectangular 
(Eulerian) form. Rectangular integration is 
effected by setting the quadrature numbers 
used for spanwise integration equal to the 
input segment lengths (loc. (100-114)). Use 
of this option is recommended for blade con- 
figurations with significant discontinuities 
in spanwise properties (e.g., counterweights, 
tip weight, step twists, etc.). 


3D 



Location 

52 

*53 

55 

56 

57 

58 

59 


Item 


(Control) 


(Control) 


(Control) 


a 


s 


L 

req 


Description 

Make nonzero (1.) to output modal integration 
constants used in the eigensolution and, to a 
limited extent, in the time-history solution. 

Make nonzero (1.) to load (vorticity) Induced 
velocity distribution either from cards or 
unit 14, as per location 999. 

Make nonzero (1.) to use the j_nt^ced velocitie s 
loaded as per location 53. 

Shaft angle-of-attack, deg. This input item 
serves a dual role: when variable inflow is 

used (either vorticity induced or Glauert 
momentum) a is used to define A . Also, 
when a major iteration is to be performed 
(nonzero location 60) and when trims on pro- 
pulsive force are deactivated (zero location 59) 
the program will trim to this input shaft angle. 

Requested value of lift to be used in major 
iteration, lb. 


PF 

req 


Requested value of propulsive force to be used 
in major iteration, lb. 


^Llft 


Tolerance on lift for major iteration, lb. 
A zero value deactivates trimming on lift. 


e 


PF 


Tolerance on propulsive force for major itera- 
tion, lb. A zero value deactivates trimming 
on propulsive force. The automatic trim 
calculation (major iteration) must trim either 
to a required propulsive force or to a required 
shaft angle-of-attack; therefore a deactiva- 
tion of trim to propulsive force automatically 
directs the trim calculation to trim to re- 
quested shaft angle, location 55. 



Location 

Item 

' Definition 

ao 

1 

' Maximum number of major iterations to be made 

in an attempt to achieve trim, A aero vajue 
will deactivate the major iteration; a nega- 
tive value will activate a stall avoidance 
calculation if, when attempting to trim, tlie 
rotor becomes stalled. (See section on genera J 
information for efficient program usage.) 

61 

, 

Built-in precone angle, deg. An identically 
zero value denotes a rotor blade articulated 
in flapping; similarly, a finite nonzero sig- 
nifies a blade nonarticulated in flapping or 
a teetering rotor. 

62 

• PM 

req 

Requested value of pitching moment for ma jot- 
iteration, Ib-f t (positive nose up) . 

63 

c 

PM 

Tolerance on pitching moment for major itera- 
tion, Ib-ft. A zero value deactivates trim- 
ming on pitching moment . r; 

64 

RM 

req 

Requested value of rolling moment for major 
Iteration, Ib-ft (positive port side up) . 

65 

^RM 

Tolerance on rolling moment for major itera- 
tion Ib-ft. A zero value deactivates trim- 
ming on rolling moment. 

66 

(Control) 

Trim deriva,tive calculation option. Input 
(0., 1.) to calculate the aerodynamic 
partial derivative matrices (using linear 
strip theory, numerically). 

67 

8 

Initial condition on articulated flapping 
or teetering angle, as appropriate, rad. 

68 

* 

6 . 

Initial condition on (nondimensional) artic- 
ulated flapping or teetering angle rate, as 
appropriate. 

69 

6 

Initial condition on articulated lead-lag 


angle, rad. (positive forward). 


52 



Location 

Item 

Definition 

70 

-k 

(6) 

Initial condition on (nondlmenslonul) 



articulated lead-lag rate. 


Built-in prelead angle, deg. An identically 
zero value denotes a rotor blade articulated 
in lead-lag; similarly a finite nonzero value 
signifies a blade nonarticulated in lead-lag. 

73 (Control) Input nonzero (1.) to harmonically analyze 

and output harmonics of flatwise stresses . 

7A (Control) Input nonzero (1.) to harmonically analyze 

and output harmonics of edgewise stresses . 

75 (Control) Input nonzero (1.) to harmonically analyze 

and output harmonics of torsional stresses . 

76 (Control) Initial condition input control. Input (1., 2.) 

to (only write, both read and write) response 
initial conditions to unit 26. Note that an 
input value of 2. is typically used only for 
restart cases involving previous runs. 

77 i|) Initial condition on rotor azimuth, deg. 

78 K . Effectivity factor of the induced velocity 

vim , . ... 

calculated using actuator disk momentum con- 
siderations in calculating inflow angle at a 
local blade section. Default value is 1., 
corresponding to conventional usage of 
momentum actuator disk inflow. This input 
quantity can be used to approximate the 
effects caused by real inflow characteristics 
as modeled by more accurate theories. For such 
usage, the effectivity would typically be in 
the range of 1.0 -*■ 1.1. 

79 (Control) Generalized Glauert (momentum derived) variable 

inflow option. A zero value deactivates usage. 
A value of 1. causes the input induced velocity 
components to be used; a value of 2. causes the 
input values to be used initially, and then 
varied in trim calculations in place of control 



5 ) 



Location 


Item 


Definition 


80-82 

83 

84-89 

84 

85 

86 
.87 

88 

89 




(Control) 


us£. 




angles; a value of 3. causes the control and 
shaft angles to be fixed and the induced vel- 
ocity components to be varied only to satisfy 
momentum equations in any requested major 
iteration. If the value is 1., the induced 
velocity components will be varied to satisfy 
momentum considerations in addition to the 
usual trim calculation. 

Initial conditions on the "momentum" induced 
velocity components comprising a Glauert- 
llke variable inflow description. Note, 
that the "vorticity" variable inflow (con- 
trolled by locations 53 and 54) and the 
momentum variable inflow can be used separately 
or simultaneously. 

Input (1.) for teetered rotor configurations. 

These input data pertain only to teetered 
rotor configurations (loc (83) = 1.) 

Distance the hub vertex point is underslung 
below the teeter hinge, in. 

Limit stop angle below which there is no 
constraint on teetering, deg. (a zero input 
is treated as infinity). 

Initial limit stop flapping spring rate defined 
when 3^ limit stop is first reached, ft-ib/rad. 

Saturation limit stop angle above which the 
limit stop constraint is defined by a second 
constant spring rate, deg. 

Saturation moment when the ^2 limit stop is 
reached, ft-lb. 


Kg Flapping spring rate of limit stop for teeter 

2 angles in excess of ft-lb/rad. (Note 

that locations 85-89 provide a general rep- 
resentation of a nonl Ine.ir limit stop spring.) 
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Location 

Item 

Description 

93 


Blade airfoil section skin friction drag 
coefficient, for use with swept airfoil, 
skewed flow option. 

94 

4v 

''se 

Critical (viscous) damping ratio used to 
approximate structural damping in edgewise 
bending modes. 

95 

(Control) 

Variable (vorticity) inflow shape option. ■ 

A value of (0.,, 1.) causes the input yorticity 
variable inflow azimuthal distributions to 
be used in a (continuous point-to-point, 
stepped pulse) manner. 

96 

(Control) 

Make 1., to invoke the radial flow swept 
airfoil option. 

97 

01 

equ 

Linear equivalent blade (nonlinear) twist 

angle defined similar to location 22, deg. 

A nonzero value for 6i is required to use 

^qu 

the input nonlinear twist arrays and is used 
to calculate that portion of the aerodynamic 
inflow velocity at the 3/4 chord position, 

U , accruing from radial flow and twist. 

P 

98 

(Control) 

Input nonzero (1.) for stress calculations 
using the mode deflection method. Zero value 
defaults to force-integration method. 

99 

(Control) 

Location used to end a case or series of 


cases. Input (+1.) to end the Loader Format 
data block for the case defined by the Loader 
data and load additional cases at the con- 
clusion of that case. Make (-1.) to end the 
Loader data and read no further cases. In 
both instances the word count, LL, (see 
beginning of this section above) should be 
(-1). Note; this entry must appear singly 
on an input card, and that card must be the 
last card for the case. 
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Location 

100-114 

115-129 

130-144 


145-159 

160-174 

175-189 

190-204 

205-219 


Item 

Ax 


M 


i 


c . 
1 



K 


t/M 



Description 

Nondimensional blade segment lengths, in order 
from root to tip, maximum of 15 values, start- 
ing from the offset location. Accuracy is 
generally improved if the last segment is small 
(< 0.03). 

2 

Mass of each blade segment, lb-sec /ft. 

Aerodynamic built-in nonlinear twist angle dis- 
tribution, deg. Since collective angle is de- 
fined at the 15 % span location, Og should have 
a zero value at 75% span. Should the structural 
twist angle distribution differ from Og the 
appropriate data must be loaded into locations 
690-704; otherwise, will be used for both 
aerodynamic and structural applications. Non- 
linear twist distributions will be used only 
if a nonzero value is input into location 97, 

®1 

•^equ. 

Blade chord at center of each blade segment 
(use for nonconstant chord blades only), root 
to tip, ft. 

Constants relating (nondimensional) flatwise 

second derivative to flatwise stress (Ec/R)„, 

F 

evaluated at center of each segment, root to 
tip, psi. 

Constants relating (nondimensional) second 

derivative to edgewise stress (Ec/R) , root 
, . E 

to tip, psi. 

Constants relating torsional moment to tor- 
sional stress, root to tip, in 

Chordwise mass radii of gyration of blade 
segments about elastic (reference) axis, 
root to tip, nondimensional with respect to 

R. 
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Location 

220-234 

235-249 

250-264 

265-279 

280-303 

280-283 

286-289 

292-294 

295-297 

298-300 

301-303 

310-333 


Item 



^lOc/4 






* 


q 


0 . 

1 


Description 

Thicknesswise mass radii of gyration of blade 
segments about axis perpendicular to chord 
line and through the reference axis, root to 
tip, nondimensibnal with respect to R. 

Area radii of gyration about elastic axis, 
root to tip, nondimensional with respect 
to R. 

Distances from elastic axis forward to airfoil 
quarter chord position, root to tip, nondimen- 
sional with respect to R. 

Distances from elastic axis forward to airfoil 
section mass centers, root to tip, nondimensional 
with respect to R. 

Initial conditions on the elastic modal re- 
sponses for the reference blade . 

Initial conditions on i'th flatwise bending 
mode deflections. 

Initial conditions on i’th flatwise bending 
mode (nondimensional) rates. 

Initial conditions on k'th edgewise bending 
mode deflections. 

Initial conditions on k'th edgewise bending 
mode (nondimensional) rates. 

Initial conditions on j ' th torsion mode 
deflections. 

Initial conditions on j ' th torsion mode 
(nondimensional) rates. 

Initial conditions on the elastic modal re- 
sponses for the second blade of a teetered 
rotor configuration, in same order as loca- 
tions 280-303. 
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Location 


Item 


Description 


355-369 

498 

499 

500 

501 

502-549 

554 

559-569 

600-614 

615-629 


A 

s 

NLDATT 


(Control) 


(Control) 



(Control) 



(I/C)j 


(1/c)^ 


Aerodynamic built-in sweep angle distribution, 
positive aft, deg. 

Segment number of the lag damper (outboard) 
attachment point. Default value is 1. 

Input (0., 1.) as the o utbo ard lag damper 
attachmen t is (inboard, outboard) of the 
feathering bearing. 

Input (1.) to use tabulated lag damper moment 
vs. angular velocity properties . 

Number of abscissa-ordinate pairs used to 
define the lag damper moment-angular velocity 
characteristics. Maximum of 24 pairs are 
allowed. 

Table of lag-^amper^charac^eristics taken in 
point pairs; (Aj^^), M(Aj^) ; U(A) = rad/sec; 
y(M) = ft-lb. (Normal sign convention is that 
A and M have the same sign.) 

In general, a nonzero value invokes the 
hysteretic (no nviscou s) structural damping 
formulation in the edgewise bending mode 
equations. Specifically, make (1., 2.) to 
use (the constant location 94 value, a span- 
wise distribution of values) in this non- 
viscous formulation. 

Distribution of edgewise nonviscous structural 
damping (= 1/2 g , where g is the usual 
structural damping coefficient). 

Section^modulii for flatwise bending, root to 
tip, in . Note that the products of (1/c) 

r 

and Kp (locations 160-174) must equal the 
bending stiffness distribution, EI/R. 

Section moduli! for edgewise bending, root 
to tip, in • See remarks above for flatwise 
bending section moduli!. 



Location 


Item 


Description 


660-674 . EB^ 


675-689 EB^ 


690-704 9g 


705-719 A0 

B 


720-734 



735-744 


745-754 GJ „ 

FB 


Torsional stiffness (to be multiplied by twist 
rate squared), as defined in reference 8, 
Ib-ft^. 

Torsion to edgewise elastic coupling stiffness 
(to be multiplied by ^wist rate), as defined 
in reference 8, Ib-ft . 

Structural built-in nonlinear twist angle 
■distribution, root to tip if different from, 
aerodynamic twist, deg. See remards above 
for aerodynamic built-in twist, location 
130-144. 

Built-in (structural) twist angle change 
per segment length distribution, root to 
tip deg. Note that this item is a direct 
statement of the built-in twist rate distri- 
bution, 0'; if all values of this distribu- 
3 

tion are xnput as zero, the twist rate distri- 
bution is computed internally using numerical 
methods from the input twist angle distribu- 
tions, locations 130-144 or 690-704, as 
appropriate. 

Distances from reference (elastic) axis 
forward to edgewise bending neutral axis, 
root to tip, nondimensional with respect 
to R. 

Flexbeam plate bending stiffness distribu- 
tion, Ib-ft. Note that all items in locations 
735-824 are required only for a redundant 
analysis of the torque-tube (flexbeam span of 
bearingless rotor (CBR) configurations) (see 
location 991) . 

Flexbeam torsion (S^. Venant) stiffness 
distribution, Ib-ft . 



Location 

Item 

Description 

755-76A 

Sb 

Flexbeam section width (chord) distribution, 
ft. 

765-774 


Torque tube flatwise bending stiffness distri- 
bution, Ib-ft . 

775-784 


Torque tube edgewise bending .stiflness distri- 
bution, Ib-ft . 

785-794 

1 

Mass of each torque tube segment whose^span 
is defined by location 100-109, lb-sec /ft. 

795-804 

(y ) 

'•^10 TT 
CG 

Distances from reference axis forwarded to 
torque tube mass centers whose masses are given 
in locations 785-794, nondimensional with 
respect to R. (Input locations provided, but 
quantity is not used.) 

805-814 


Section mod^lii for torque tube flatwise 
bending, in . (Input locations provided, 
but quantity not used.) 

815-824 

(I/c)ett 

Section mod^lii for torque tube edgewise 
bending, in . (Input locations provided, 
but quantity not used.) 

850-881 

[T] 

Trim partial derivative matrix, (obtained 
as output from a previous trim calculation. 
Nonzero values will suppress internal re- 
calculation of [T] . 

900 

^pend 

Number of blade appended pendular vibration 
absorbers activated, maximum of 2. 

900-910 

- 

Data for pendular absorber No. 1. 

901 

’^PAH 

Radial location of pendular absorber hinge 
attachment, in. 

902 


Location of pendular hinge point in direction 
normal to the section chord, (+ upward), in. 

903 


Pendular absorber arm length (i.e., distance 
from hinge point to mass center), in. 
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Location 


Item 


904 

905 


M, 


906 


K 


907 

908 

909 


8 


P 

( 1 ) 


910 


*(I) 

P 


911-920 


921-922. 


PA 


923-924 


925-928 


PA 


925-926 


P Pi 


Description 

2 

Mass of pendular absorber, lb-sec /ft. 

Rotary inertia of pendular absorber about . 
the hinge point, lb-sec -ft. 

Discrete rotary spring restraining the 
pendular absorber about the hinge point, 
ft-lb/rad. 

Rotary damping about hinge point, Ib-sec-ft. 

Frictional moment about the hinge point, Ib-ft 

Initial condition on pendular absorber flap 
angle, rad. 


Initial condition on pendular absorber 
(nondimensional) flap rate, rad. 

Repeat of locations 901-910, but for the 
second pendular absorber . 

Integration frequencies for the pendular 
absorbers 1 and 2, respectively. These 
frequencies should be close to the respec- 
tive natural frequencies of the absorbers. 
Default values are calculated from the input 
mechanical properties, nondimensional. 

Built-in (+, nose up) pitch angles, relative 
to section chord line for pendular absorbers 
1 and 2, respectively, deg. 

Initial conditions for the pendular absorbers 
of the second blade of a teetered rotor. 

Initial conditions for pendular absorber No. 1 
(same as locations 909 and 910, but for second 
blade), rad. 



Location 


927-928 

930-939 

930 

931 

932 

933 

934-936 

939 

940-945 


Item 


6 


( 2 ) *( 2 ) 
P2 ’^P2 


m 


Aa 


Aa 


NSTAVB 


ISTAVB 


(Control) 


Description 

Initial conditions for pendular absorber No. 2 
(same as locations 919 and 920, but for second 
blade), rad. 

Input data for controlling the calculations of 
the hub excitation vector and the impedance 
matrix . 

Multiple of number of blades to determine the 
harmonic of hiib motion and resulting loads 
and impedance. Thus, resulting frequency, 
w, would be u) = mbf2. 

Perturbational amplitude used for incrementing 
each of the components of linear hub accelera- 
tion (x,y and 2) to calculate the impedance 
matrix numerically, ft/sec . 

Perturbational amplitude used for incrementing 
each of the components of angular hub accelera- 
tion (6„,e ,6_) to calculate the impedance 
matrix numerically, rad/sec^. 

Number of blade radial stations for which 
harmonics of vibratory stresses are to be 
included in augmented impedance matrix, 
maximum of 3. 

The (NASTAVB) radial station segment numbers 
for utilizing vibratory stresses in the aug- 
mented impedance matrix. 

Input (1.) to include vib rator y p ush-rod loads 
in the augmented impedance matrix. 

Higher harmonic control pitch angle amplitudes 
wherein the frequencies are at (b-1) , (b) and 
(b+1) times the rotor rotational speed and b, 
the number of blades, is location 5, and where 
a negative Fourier Series representation is 
assumed. 
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Location 

Item 

Description 

940-941 

^(b-1) ’®(b-l) 
s s 

Cosine and sine components of (b-1) per rev 
higher harmonic control pitch angles, respec- 
tively, deg. 

942-943 


Cosine and sine components of (b) per rev 
higher harmonic control pitch angles, respec- 
tively, deg. 

944-945 

^(b+1) ’®(b+l) 
s s 

Cosine and sine component.s of (b+1) per rev 
higher harmonic control pitch angles, respec- 
tively, deg. 

946 

■ A0 

HHC 

Perturbational amplitude used for incrementing 
the components of higher harmonic control 
(locations 940-945) to calculate the augmented 
impedance matrix numerically, deg. A zero 
value will retain the input higher harmonics 
of pitch control, but will suppress this por- 
tion of the impedance matrix calculation. 

957 

h 

Height of positively thrusting rotor from 
ground or wind tunnel floor for purposes of 
evaluating Heyson corrections to rotor angle 
of attack (see Ref. 5), ft. 

958 

»WT 

Height of wind tunnel test section for purposes 
of evaluating Heyson corrections, ft. Note 
that zero values for the wind tunnel test 
section dimensions implies that ground effect 
corrections, rather than wind tunnel wall 
corrections, are to be made. 

959 

”WT 

Width of wind tunnel test section for purposes 
of evaluating Heyson corrections, ft. 

975 

CASE 

Case number. 

976 


Built-in elevation angle of push-rod attach- 
ment point, measurable when 6 = 0, deg. 

977 

NDEL3I 

Segment number of inboard attachment point 


of pitch input structural members (pitch horn, 
feathering cuff or torque lutie) to blade .spar. 



Location 

978 


979 

980 

981 

982 

983 

984 


t 


Item 


NDEL30 


r 

PR 



NSPSTR 


GJ 


FB 


Description 

Segment number of outboard attachment point 
of pitch input structural member to blade 
spar. Note that if this structural element 
is attached to the blade spar at a slngJe 
point (cantilevered configuration) NDEL3J 
and NDEL30 must both be Input with the same 
appropriate value. For bearingless rotor 
(CBR) applications NDEL30 serves the additional 
function of defining the. outer limit of the 
flexbeam; this limit is taken as the inner 
boundary of the NDEL31'th segment. 

Radial location (blade station) of push-rod 
to pitch input structural member' attachment 
point, in. 

Location, .forward of feathering axis, of push- 
rod to pitch input structural member attach- 
ment point, in. Note that input items 977- 
980 together with the input modal data provide 
the geometric data from which pitch-flap, 
pitch-flat and pitch-edge coupling are calcu- 
lated. These calculations are bypassed if 
either of locations 978 or 980 are input 
identically zero . 

Segment number of innermost end of the feather- 
ing flexure; defaults to 1. 

Built-in twist angle of the flexbeam, positive 
leading edge up at outboard end of flexbeam, 
deg. Note that •^Bpy defined relative to the 
inputted twist angle distribution, which in 
turn is defined for zero collective angle. 

Thus, Gg represents the built-in collective 

X* D 

angle for the blade. 

(Inverse) averaged torsional stiffness of the 
flexbeam Ib-ft . 


k Equivalent or critical torsion stress coefficient 

FB for flexbeam; i.e., stress per unit average 

torsion strain as defined by total flexbeam 
twist, in-psi. 
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Location Item 


Description 


985 (Control) 


986 (Control) 


987 K„ (K „) 

= 2 ^ Sp 

988 Kg (Kgp) 

Sys S£ 

989 (Control) 

990 NTQTBl 


Make nonzero (1.) to include the "wobble " 
mode for cantilever torque tube configured 
bearingless rotors. 

When finite snubber stiffnesses . are input 
(locations 987 and 988), the stiffness 
alignments (principal axes) are assumed 
to be (vertically and horizontally, in 
torque tube flatwise and edgewise directions) 
as location 986 is (0., 1.). 

Snubber vertical (torque tube flatwise) 
spring rate, Ib/ft. 

Snubber horizontal (torque tube edgewise) 
spring rate, Ib/ft. 

When the torque tube flexbeam redundant 
analysis is involved (location 991) the 
torque tube is assumed to be in ( tension , 
compression ) as location (989) is (0., i.). 

Segment number of innermost end of torque 
tube; defaults to 1. 


991 (Control) 


992 (Control) 


993 

FB 


Input nonzero (1.) to activate the redundant 
analysis option. When the redundant analysis 
is activated, the flexbeam and torque tube 
stiffness and mass inputs, locations (735) 
through (794), must be appropriately nonzero. 

Used to weight the effectiveness of nonlinear 
AEI bending excitation of torsion over the 
flexbeam span; the effectiveness is taken to 
be (100%, 0%, 100% (x) pseudo-torsion mode 
shape) as LOCATION (992) is (0., 1., 3.)). 

Flexbeam aspect ratio parameter (=- /l.5(l-v)) 
to determine plate-like deflection for torsion 
"pseudo-mode". A zero value gives a rod-like 
deflection shape. 
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Location 


994 


Item 


995 




996 (Control) 


997 (Control) 


998 


Ai|) 


F389 


999 (Control) 


Description 

Spring rate of rotational spring connecting 
outboard end of torque tube to the blade 
spar about chordwise axis, ft-lb/rad. A 
zero value implies zero stiffness for 
pinned-pinned torque tube configurations 
and infinite stiffness for snubbed config- 
urations . 

Spring rate of rotational spring connecting 
the outboard end of the torque tube to the 
blade spar about an axis perpendicular to 
the chord line, ft-lb/rad (see above). 

Analytic (static) airfoil option . Make 
nonzero (1.) to use the built-in analytic 
approximation to the static NACA 0012 air- 
foil data. 

Unsteady airfoil data option . Make 
nonzero (> 1.) to use the synthesized 
unsteady airfoil data technique described 
in Reference 7. This can be used with 
either the explicity input airfoil data 
or the built-in analytic NACA 0012 static 
airfoil data. For nonzero values greater 
than unity, the "cutoff Mach number" (Mach 
number above which the unsteady data is 
disregarded in favor of quusi-statlc data) 
is taken to be the input value minus one 
(1.) rather than the built-in default value 
of 0.6. However, only nonstandard cutoff 
Mach numbers between 0.1 and 0.95 are accepted. 

Aximuth increment used to generate punched 
card image output which is then used as input 
data to UTRC Program F389 for computing 
variable inflow. This option is bypassed 
with a (0.) input value. 

Input (0., 1.) for data transferrals with the 
F389 variable inflow program to be made using 
(punched cards, data files). If data files 
are used, G400PA reads from unit 14, and 
writes to units 16, 23 and 13 . 
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Location 


Item 


Description 


1000 

1001 

1002-1050 

1051 

1052-1100 

1101 


(Control) Input nonzero (1.) to activate use of the 

tabulated time-histories of incremental 
control angles. 

Number of abscissa-ordinate point pairs 
used to define time-history of A0 ^^(t) ; 
calculation of this time-history is bypassed 
with a (0.) value. 

A0 Table of AO abscissa-ordinate pairs; 

U(A0 <= <3eg.; U(t) = sec. . 

NAAj^ Number of abscissa-ordinate point pairs 

used to define time-history of AAi (t) ; 

g 

zero value bypasses calculation. 

Table of AA abscissa-ordinate point pairs; 

U(A-| ) = deg? U(t) “ sec. 

■^s 

NABj^ Number of abscissa-ordinate point pairs 

® used to define time-history of AB^^ (t) ; 

zero value bypasses calculation. ® 

ABj^ Table of ABj^ abscissa-ordinate point pairs; 

® U(Bi ) = degf U(t) = sec. 

s 




1102-1150 



Last Card Image for Block II (Loader) Data : 


f -1 99 ± 1. 


(12, 14, F12.0) 


other cases follow 

preceedlng data defines last case 


111. Blade Mode Shape Data 

Included in this data are the radial distributions of the blade 
(uncoupled) flatwise, edgewise and torsion normal mode shapes and their 
derivatives. These quantities must generally be provided from an external 
source such as United Technologies Corporation Program E159 or an equivalent, 
in the following card image format: 


^ NFM 

NEM NTM 

NSEG 

(414) 

subsequent 

cards : 



^ F(I) 

F(I+1) 

F(I+2) F(I+3) F(l+4) 

(F18.0, 4F12.0) 


where: NFM, NEM, and NTM are, respectively, the numbers of flatwise bending, 

edgewise bending and torsion normal modes whose mode shapes and derivatives 
are to be input. NSEG is the number of blade spanwise stations for which 
the input modal data are defined. F(i) are the modal functions listed below 
(defined at the i'th spanwise stations). Five entries per card are made 
for each F function input for NSEG total entries. The modal functions must 
be loaded in the following order: 
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(NSEG values, root to tip) 



( " 


It 


It 


" ) 


(NFM flatwise modes) 




(NEM edgewise modes) 




(NTM torsion modes) 



IV. Var j ;jbi(i ■ llurmonJ c Inflow 

If location 53 of the Loader block of operational data is nonzero, the 
following block of variable harmonic inflow is loaded in: 



NHARML is the number of harmonics of inflow to be loaded, and XLAMO(i) , 
XLAMC(i,n) and XLAMS(i,n) are, respectively, the zeroth, n'th cosine and 
n'th sine components of harmonic inflow at the i’th radial station, wherein 
a positive Fourier series is assumed. The inflow is defined positive up and 
has the units of ft/sec. If n > 5 card image #3 is repeated until n cosine 
coefficients are input. The series of card images, #1 through 1 ) 3 + are then 
repeated for the remainder of the NSEG segments. 

V. Multiple Case Runs 

The above described data set-up defines the correct ordering of required 
data blocks for a general case. When multiple cases are run (while remaining 
within the G400PA branch of the base program) the second and subsequent cases 
utilize most of the data input for the first case. The following rules apply 
to the running of multiple cases: 
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1. Airfoil data is loaded only. for the first case; all subsequent cases 
within the run use the same tabular data, if analytic data is not used, 

2. Only those items within the operational (Loader) data whicrh are to be 
changed from case to case need to be input. 

3. Item 99 of the operational data controls the running of subsequent 
cases; a (+1.) value causes a subsequent case to be loaded whereas a (-1.) 
value terminates the computer run after the current case. 

h. Unless otherwise specified (by a +1. value for operational data item 
19) the input modal array data block is used for all cases within the run. 

5. Similarly, unless otherwise specified (by a +1. value for operational 
data item 53) the input harmonic variable inflow data block is used for all 
cases within the run. 

6. Operational data items 19 and 53 discussed above are both automatically 
set to zero at the conclusion of the data input for every case. 

7. Terminal conditions on the blade azimuth angle, item 77, and on the 
degrees-of-f reedom, items 67-70, 280-303, 310-333, 909, 910, 919, 920, and 
925-928, for any case are carried over as initial conditions on these quan- 
tities for the subsequent case. Thus, for some applications, e.g., investiga- 
tions of unstable responses, it would be appropriate to reinitialize these 
items on the subsequent cases. 

When solution flow leaves the G400PA branch of the base program, the 
ability . to carry over terminal conditions (as initial conditions for a 
subsequent case) and/or any other quantities associated with trim is lost. 
However, a need still exists for preserving these initial conditions for 
subsequent reentries to the G400PA branch of the program. As per loader 
locations (76) (described above) these initial conditions are written to 
and read from file 26. 



General Information to Facilitate 
Operation of Program and Improve Efficiency 

Aside from considerations of the actual aeroelastic parameters describing 
the blade configuration, which are covered in the above sections, additional 
attention should be paid to the mechanics of obtaining efficient numerical 
solutions of the dynamic equations. In this regard, there arise two basic 
areas of concern wherein this section should be of assistance. The first of 
these, areas is the proper selection of parameters for efficient temporal 
numerical integration of the dynamic equations (flapping or minor Iterations) 
and the second is the proper selection of parameters for effecting a satis- 
factory rotor trim (major iteration). The following subsections provide 
information for making proper parameter selection in eacli of tiiese areas. 

Temporal Numerical Integration 

As discussed in Reference 1 temporal integration of the higher differen- 
tiated response variables to obtain the lower ones is achieved in the G400PA 
program using a variant of the Adams integration algorithm. The selected 
algorithm is defined by means of the azimuthal integration step size, A\j;, 
and the integration frequency, oi. 

The integration. step size should be an integral divisor of 360; a proper 
choice depends on the maximum coupled frequency inherent in the various aero- 
elastic responses. A reasonable upper limit for AiJ; is 30 divided by the 
maximum such frequency in per rev. Values of A>Jj greater than this upper limit 
will compromise the integration accuracy and, for sufficiently large values, 
will cause the computed responses to develop "numerical" instabilities. As 
a corollary, a check on any response which is predicted to be unstable by 
the analysis, is to rerun the case with a reduced integration step size to 
test for the possibility of the unstable response being merely a numerical 
instability. 

For each response degree-of-freedom a different integration frequency, 

( 0 , is used in the integration algoithm; this frequency is, for. each of the 
elastic modes, the respective input natural frequencies (locations 31-34, 
37-39, and 40-42). The usual pendular frequencies of a rotating, hinged, 
rigid beam are taken as the integration frequencies of the flapping and 
lead-lag degrees-of-f reedom. In addition to defining modal stiffnesses and 
integration frequencies, the input frequencies serve yet another purpose. 

As noted above, the proper value of integration step size, Ai|;, varies in- 
versely with the maximum modal frequency. Thus, run times (caused by reduced 
step size) will significantly increase as any one modal frequency increases. 
Since any degree-of-freedom exhibiting a large natural frequency tends to 
respond quasi-statically , i.e., as if the acceleration ( q ) term were 
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negligible, a reasonable approximation the response calculation is to 
avoid the numerical integration of the q term entirely and treat the 
response quasi-statically . This option can be invoked for any such high- 
frequency mode by input of a negative frequency; a negative sign will not 
affect the proper usage of the frequency in the calculation of the 
dynamic equations. Note that this optional response calculation can be 
invoked singly or in combination for any of the el.'istic modal responses 
(negative values in any of locations 31-34, 37-42, 921 and 922). 

Hub Force and Moment Trim 

Operation of the trim or major iteration feature of the program is 
controlled by input locations 55 through 60, 62 through 66 and, for some 
applications, 79. The main control for the major iteration is location 
60, the number of major iterations, ^ zero value causes the major 

iteration feature to be completely deactivated. On the basis of past 
usage, a reasonable range for this input appears to be from 5 to 10, 
depending on the tightness of the convergence tolerances selected and 
the "goodness" of the initial guesses on the control parameters. Conver- 
gence of the major iteration is adversely affected by any lack of con- 
vergence of the responses, to periodicity within each trim iteration and 
by incursion of the rotor into a significantly stalled flight regime. 

Note that a negative value of N will activate a stall avoidance calcula- 
tion wherein the controls will be perturbed to reachieve unstalled flight 
rather than to achieve the prescribed trim condition. Should a major 
iteration fail to converge within any one run the last used control 
angles and Initial conditions on response variables are generally avail- 
able in output card image form and/or partially in the output printed 
records of each major Iteration for use in subsequent major iterations. 

An additional capability available for achieving a trimmed configura- 
tion is the option to calculate the aerodynamic partial derivative matrix 
either approximately using simple linearized strip theory or using a 
straightforward numerical differentiation. This option is controlled by 
loader location 66. It should be noted that the accuracy of the numerical 
differentiation is dependent upon how well the responses have converged to 
periodicity, which in turn is controlled by the flapping tolerance, loc. 11. 
When numerical derivatives are calculated, they are output on card images 
for subsequent input, using loader locations 850-881. 

The trim iteration is operationally flexible as to what hub loads it 
will drive to requested values. Generally, the various reciuested hub 
loads, lift, propulsive force, pitching and rolling moment are activated 
in turn by specifying nonzero values for each of their respective tolerances. 
Specifically, the following table describes the optional combinations of 
hub loads and rotor shaft angle available with the G400PA trim capability. 



TABLE II 


SUMMARY OF OPTIONAL BASIC TRIM COMBINATIONS 


Option 

Lift 

Prop . 
Force 

Pitch. 

Moment 

Roll 

Moment 

®.75R 

^1 

s 

S 

a 

s 

1 

S 

U 

U 

U 

V 

F 

F 

S,(V) 

2 

S 

S 

U 

U 

V 

F 

F 

V 

3 

U 

u 

s 

s 

F 

V 

V 

•S,(V) 

4 

s 

u 

s 

s 

V 

V 

V 

S,(V) 

5 

s 

s 

s 

s 

V , 

V 

V 

V 

6 

u 

u 

s 

U 

F 

V 


S,(V) 

7 

s 

u 

s 

u 

V 

V 


S,(V) 

8 

s 

s 

s 

u 



V 

V 


V 




where: F: Control parameter kept fixed 

V: Control parameter varied 

S: Trim parameter specified and trimmed to 

U: Trim parameter unspecified and ignored 
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Note that for each option available, the individual trim parameter is 
activated by input of a positive, nonzero value for the respective parameter 
tolerance. For each option the trim calculation trims either to the specITied 
propulsive force or the input sliaft angle, a (location 55) as tlie propulsive 
force tolerance, is finite or zero, respect Lvei y . A reasonabh: eliolce 

of tolerance values for the hub force and moment loads can be obtained I rom 
considerations of the helicopters inertia properties and acceptable devia- 
tions of the load factor from unity, and of the pitching and rolling accera- 
tions from zero. 

The trim calculation can also be used in conjunction with the Glauert 
variable inflow. Location 79 controls the use of this simplified form of 
variable inflow. A variety of trim calculation operations are achieved 
when this control location is nonzero. If this control option is 1., the 
Glauert variable inflow induced velocity components- v^, and v^ 

(locations 80-82) would be varied in addition to those parameters s^own in 
Table II to effect momentum balances in thrust, pitching moment, and rolling 
moment as well as the specified trim. If the control option is 2., only the 
velocity components would be varied (instead of the control angles, 0 
A , and but in an analogous manner. If the control option is 3., 

tnl trim cafculation would again only vary the velocity components, but to 
achieve momentum balance only, with an specif ied • trim. 
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PROGRAM OUTPUT DESCRIPTION 


The complete printed output generated by the G400PA program can be classified 
into the following four major categories: 

I-. Listing of Input Data 

II. Parameters Calculated from the Input Data 

III. Results of Time-History Solution 

IV. Impedance Matrices 

This section describes the pertinent output pages associated with each of these 
categories. It should be noted that while output will always be generated for 
the first three categories, output for the remaining category depends upon the 
optional calculation of these impedance matrices. The subsections which follow 
describe, in turn, the details of each of these four categories. 


Listing of Input Data 

Output in this category includes a listing of the static airfoil data (if any 
is inputted), a descriptive listing of the "inertia, elastic, geometric and other 
operational" (Loader Format) data, and a listing of the harmonic variable Inflow: 
blocks I, II, and IV, respectively, of the input data. If static airfoil data is 
inputted, then a listing of this data will be outputted for c^, cj and c^ m each 
with the format shown in Sample Page 1 where each column represents data at one 
Mach number. Within each column the first line gives the number of angle-of-attack/ 
aerodynamic coefficient pairs defining the functionality; the second line is the 
Mach number, and the ensuing line pairs are the angle-of-attack/aerodynamic coeffi- 
cient pairs, where the angles-of-attack are in degrees. This output closely follows 
the input format described in the previous section. 

A description of the Loader Format data output is omitted herein since this 
output merely duplicates the description already given in even greater detail in 
the previous section. In Sample Page 2 is showTi the listings of the inputted 
harmonic variable inflow, both by harmonic components and by azimuthal variation 
for each of the (maximum of) fifteen blade segments. As with the input format, 
the inflow is defined positive up, has the units of ft/sec and a conventioanl 
positive Fourier series representation is used. 


Parameters Calculated from the Input Data 

Sample Pages 3 and 4 list a variety of intermediate parameter calculations which, 
for the most part are directly applicable only to the composite bearingless rotor. 
Sample Page 3 presents the matrices and vectors used to define, respectively, the 
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elastic bending characteristics of the torque tube and the inboard (snubber) 
end vertical deflection of the torque tube. For each, a quadratic variation 
with pitch angle is assumed so that the total bending stiffness matrix, FDEFL, 
and inboard deflection vector, Z5I, are formed in the indicated manner wherein 
the total blade pitch angle at the 75 percent span, TTHO, is taken in radians. 
Furthermore, all the results output for these two quantities are nondimen- 
sional (see the superscripts section of the List of Symbols for the proper 
nondimensionalizatlon) . The matrix FDEFL duplicates matrix S, whose rows 
and columns are defined by equation (108) of Reference 1; vector Z51 shares 
the same columnar dependence with matrix FDEFL. Finally, at the bottom of 
Sample Page 3 are evaluations of FDEFL and Z51, using the quadratic represen- 
tations, at a pitch angle equal to the inputted collective angle, It 

should be noted that the functionalizations of these quantities to quadratic- 
form is accomplished by matching exactly the quantities calculated with pitch 
angles equal to the collective angles. 

The first group of output parameters given in Sample Page A consists of 
the (nondimensional) coefficients defining the nonlinear torsional stiffness 
characteristics of the flexbeam; they correspond directly to the coefficients 
defined in Reference 1. Again the standard nondimensionalizatlon is used. 

The second group of output parameters on the sample page consists of the quad- 
ratic functionality coefficients for the pitch-flat and pitch-edge couplings 
and the angular motions of the torque tube relative to the flexbeam at the 
juncture for each of the selected flatwise and edgewise modes. The quadratic 
functionality on pitch angle duplicates that used for FDEFL and Z51 described 
in the above paragraph. Note that RELMW and RELMV are calculated only on the- 
basis of an infinitely rigid torque tube and, hence, are omitted wlien the re- 
dundant analysis option is invoked. The third group of output parameters are 
the effective torsional springs ^ 0 ^, ^>nd as defined in Reference 1, 

and the flexbeam stiffness, GJ and Tkg, respectively. The units of the springs 
are Ib-ft/rad and those of the stiffnesses are Ib-ft^. 

In Sample Page 5 are shown typical modal information for the inputted 
flatwise and edgewise bending modes. For each such mode the (nondimensional) 
modal frequency, pitch-flat (or pitch-edge) coupling and Inputted mode shape 
and spanwise derivatives are listed. In addition the listing presents the 
derived Incremental deflection vectors which account for blade twist (see 
equations (7-10), of Reference 1. Within a flatwise modal information group, 
the DVB and DVE arrays correspond to those first order Av spanwise functions 
due to built-in twist and torsional modal twist, respectively. The DWWBB, DWWBC 
and DWWCC arrays correspond to those second order AW functions due to the com- 
binations of built-in twist with itself, built-in twist with pseudo-torsion 
mode (control) twist, and pseudo-torsion mode twist with itself, respectively. 
The various arrays, DV2BP, DV2KP, DWW2BBP, DW2BCP and DWW2CCP are the first 
spanwise derivatives of the second compom-i.ts (those with superscript " 2 ") of 
the above discussed arrays, DVB, DVE, DWV.'bH, DWWBC, and DWWCC, respectively. 
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The UWE nonlinear deflection arrays correspond to the bracketed integral 
function defined in Reference 1. Within the edgewise modal Information group, 
the various arrays, DWB, DWE, DVBB, DWBC and DVVCC, etc., correspond to 
similarly defined spanwise functions involving twist and the edgewise modal 
deflection and spanwise derivative arrays. 

The output parameters presented in Sample Page 6 consist of the inputted 
torsion modal arrays together with the derived pseudo-torsion mode shape (as 
defined in Reference 1) and spanwise derivative. Tlie output parameter grouji 
presented in Sample Page 7 consists of the spanwise distributions of various 
pertinent aeroelastic quantities. The X and XCEN arrays are the nondlmcnslonnl 
distances of the centers of the segments from the offset and rotor axis, respec- 
tively. The units of the CHORD array are feet. The angle of attack descrip- 
tors THETA-AERO, PHI and ALPHA are, respectively, the geometric aerodynamic 
pitch angle, the inflow angle, and the resulting section angle-of-attack, all 
in degrees. These angles are calculated using the input initial conditions 
on azimuth angle and on the response variable deflections and velocities. The 
resulting Mach number and aerodynamic coefficients are used to define the per- 
turbational airloads used in the eigensolution. The quantity KAPPA/U is the 
spanwise variation in aerodynamic moment damping coefficient which when multi- 
plied by the local pitch rate approximates' the potential flow unsteady pitching 
moment coefficient. The quantity (Y10C/4)/C is the spanwise distribution of 
quarter chord offset from the reference axis nondimensionalized by chord. The 
dynamic and structural quantities in the third group of this output page in- 
cludes the QUAD array which constitutes the Integration weighting numbers for 
spanwise integration. The THETA-STR array is the pitch angle distribution of 
the structural principle axes and has the units of degrees and, in general, 
differs from the aerodynamic pitch angle distribution. The two arrays, 

TWIST-BLT and TWIST-TOT, are the nondimensional structural twist rate distri- 
butions of the built-in twist and, the total twist (including elastic response 
and control Inputs), respectively; these arrays have the units of radians. The 
quantities TENSB, EIYB, EIZB, and MASSB are, respectively, the blade tension, 
flatwise bending stiffness, edgewise bending stiffness and mass distributions, 
all nondimensionalized in the standard sense. The (Y10NA)/C and (Y10CG)/C 
arrays are, respectively, the edgewise bending neutral axis and mass center 
offset distributions nondimensionalized by chord. 


Results of Time History Solutions 

Sample Page 8 presents a variety of response and load quantities defining 
the transient , aeroelastic responses. The first row of parameters following 
the page title presents, for the time-history solution, a listing of the param- 
eters defining the flight condition. These parameters consist of the various 
control angles (in degrees), the inflow and advance ratios, and the nondimen- 
sionalized values of the "momentum" induced velocity components. The remainder 
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of Sample Page 8 eomprlses the typical azimuthal listing; tills listing Is 
outputted for every azimuth angle which is a multiple ol the print azimuth 
Increment, input item no. 23. 

The first of the four groups of result quantities on this sample page 
lists the spanwise distributions of the pertinent aerodynamic quantities. The 
inflow and total section angles of attack PHI and ALPHA, respectively, are in 
degrees. The MACH NO., CL, CD, and CM are self-explanatory and nondimensional . 
The airload distributions in the Z5 and Y5 directions, SAZ5 and SAYS, respec- 
tively, have the units of Ib/in. , the aerodynamic pitching moment distribution, 
MAX5, has the units of in.-lb/in. The quantities SDZ5, SDY5 and MDX5 are 
"semi-dynamic" load distributions. These distributions are dimensionally 
similar to those above described resulting from aerodynamics, but instead arise 
from all the dynamic effects except thc^^ li^j^lving second time differentia- 
tions of the response variables (i.e., q„ . . ^Ivi, » etc.). The quantity MEX9 is 
the elastic torsion moment distribution which consists of those torsion 
couplings arising from AEI, the tension-neutral axis offset and other twist 
related elastic effects; it too, has the units of in.-lb/in. 

The second group of result quantities consists of the instantaneous 
generalized excitations, XI, and the generalized accelerations, velocities, 
and displacements (Q**, Q*, and Q, respectively), for each of the modal re- 
sponse variables selected, all appropriately nondimensionalized . The general- 
ized excitations are defined to be the elements of the right hand sides of the 
dynamic equations. 

The third group of result quantities are comprised mainly of blade deflec- 
tion and stress distributions. The vertical and inplane deflections are those 
in the Z5 and Y5 directions, respectively. All stress quantities have the 
units of Ib/in. , whereas the torsion moment has the units of lb-in. It should 
be noted that, over the f lexbeam-torque tube span, the flatwise and edgewise 
stresses outputted are those only for the flexbeam whereas the torsion moments 
and stresses outputted are those only for the torque tube. The last two 
columns arc the distributions of the A and B parameters needed to define the 
unsteady airloads. The fourth group of result quantities at the bottom of the 
sample page consists of miscellaneous deflection, load and stress results for 
the flexbeam and push-rod. 

The line titled SPAR/FLEXURE PARAMETERS presents similar stresses and 
torsion moments for the flexbeam at the indicated spanwise location; note, 
however, that, the output torsion moment is for the flexbeam immediately 
inboard of the juncture. Of the remaining four quantities, PUSH-ROD 
(RELATIVE) DEFL. (in.), TORQUE TUBE ROOT DEFL. (in.), and TORQUE TUBE ROOT 
SHEAR (lb) pertain only to cantilevered torque tube configurations wherein the 
"wobble mode" option is Invoked (input location 985). All three quantities 
are defined in the positive Z5 direction. The quantity I’USII-ROD LOAD (Ib) 
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is the upward (+Z5) directed load the push-rod exerts on the pitch horn/ 
push-rod attachment point; this quantity is calculated for all blade typos. 

For composite bearingless rotors, the push-rod load accounts for the total 
blade torsion moment at the root less that torsion moment resisted by the 
flexbeam. 

After the time-history solution has either converged to periodicity or 
run to maximum flapping trials (input location 10) various integrated loads 
are calculated for one final blade revolution to form the aerodynamic per- 
formance and stress results depicted in Sample Page 9. For each of eight (8) 
performance quantities results are presented in nondimensional coefficient 
form, in nondimensional form divided by solidity, and in actual dimensional 
form. Note that ten (10) dimensional quantities are listed and the units are 
lb for forces and Ib-ft for moments, as appropriate. The quantity EQU. DRAG 
(lb) represents the combined power expended by the rotor due to rotor rota- 
tion (torque) and translation (drag) divided by flight speed. 

The next line duplicates the parameters defining the flight condition 
and Includes four (4) additional quantities which depend on the integrated 
performance for evaluation. At the beginning of the time-history calculation 
it is not known which part of the inflow ratio being used is due to ram 
effects and which due to momentum induced effects. Once the integrated rotor 
thrust is calculated, however, the induced portion of the inflow can then be 
calculated using the simple usual momentum formula derived for flight in an 
infinite continuum. The complementary portion of the inflow represents the 
ram effect from which the shaft angle-of-attack ALPHA S, in degrees, can be 
calculated. The quantity VEL ACT. is the. actual forward flight velocity, in 
knots, consistent with the advance ratio used and the shaft angle of attack. 
For finite forward flight speeds EQU. L/D is the lift divided by the equiva- 
lent drag; for hovering cases this quantity is the figure of merit. PAR. 

AREA, the rotor parasite (drag) area, in square feet, is the rotor drag 
divided by dynamic pressure. The line titled CORRECTIONS DUE TO WIND TUNNEL 
WALL INTERFERENCE: consists of recalculations of those quantities which de- 

pend on the induced portion of the inflow wherein the Induced inflow is cal- 
culated using the formulae derived for flight in wind tunnels or in ground 
effect. The remainder of Sample Page 9 consists of reductions of the various 
stresses (Ib/ln.^) and the push-rod load (lb) to median and 1/2 peak-to-peak 
values over the final rotor revolution. 

Should major (trim) iterations be used (see descriptions of input items 
60 and 62 through 65) two similar types of page output will be generated. 

The trim calculation first calculates a trim matrix of partial derivatives. 

For this calculation the blade pitch control angles (or, optionally, com- 
ponents of Glauert inflow) and shaft angle are systematically Incremented to 
form the numerical partial derivatives. During the calculation of these 
partial derivatives, the time history solution is run separately five times to 
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establish the base values and the four perturbations. An output page is 
generated for each of these five runs. Each of these pages has the combined 
form of the first three lines of Sample Page 8 (giving header Information) and 
the aerodynamic portion of Sample Page 9 (exclusive of the stress results). 

Once the partial derivative matrix has been calculated from the perturba- 
tional calculations, the second type of trim related page output is generated, 
and is shown in Sample Page 10. This page, which is repeated for each trim 
iteration, consists of the same output as is generated while incrementing the 
control quantities (as discussed above) together with additional specific trim 
iterative information. The first line of this additional information consists 
of the iteration number, the zeroth, first cosine and first and second sine 
harmonics of first flatwise mode response, in radians, and an estimate of an 
effective angle-of-attack on the retreating blade side {\p= 270°), in degrees. 
The nonzero elements of the TRIM PARTIAL DERIVATIVE MATRIX give, for each row, 
the partial derivaties of the four trim quantities (C^, Cpp, Cpj^, Cpj.{, respec- 
tively) with respect to the four control quantities being used (0 ■75, A-j^g, 

and (sin Qg), or v^, v^g, and (sin cig)), for each respective column. 

The elements of this matrix are formulated using approximate linear strip 
theory (with reversed flow effects) or numerical differentiation and are cal- 
culated for either set of control quantities. The ERROR VECTOR consists of 
the differences between the four requested trim quantities and those achieved 
in the preceeding time-history. The two lines depicted give the error vector 
in dimensional (lb and Ib-ft) and nondimensional forms, respectively. The 
TRIM VARIABLE CORRECTION VECTOR consists of those changes to the control quan- 
tities which should null the above described error vector. The correction 
vector is obtained from the premultiplication of the inverse of the G matrix 
with the error vector, but the corrections are scaled, if necessary, to prevent 
control changes of more than 2 degrees within any one iteration. The control 
parameters whose increments are depicted in this output page are, in respec- 
tive order: 6 . 75 » ^ls> ®ls» (sin cig). A, Cj, v^, and the first 

three have units of degrees and the remainder are dimensionless or nondimen- 
sionalized. 
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PART II. TIME HISTORY 
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Once the time-history solution has converged to periodicity and all ma^’t 
iterations have been completed, the program optionally performs harmonic, anal- 
yses of the azimuthal variations of various response quantities (see descrip- 
tions of input items 48 and 73 through 75). The outputs of these harmonic 
analyses are depicted in Sample Pages 11 through 13. In each of these sample 
pages the harmonic information for each response variable is contained in the 
appropriate horizontal band of five rows. The harmonics are listed by columns 
up to a maximum of 10 harmonics. All harmonic analysis output depicted on 
these sample pages assume a negative harmonic content form in keeping with the 
(negative) harmonic form conventionally assumed for the blade pitch control 
and rigid flapping angles. For each harmonic of response variable five 
quantities are outputted; these quantities are, respectively, the cosine and 
sine components, the equivalent amplitude and phase angle, and lastly, the 
amplitude of the harmonic relative to all the other harmonic amplitudes out- 
putted. Sample Page 11 depicts the harmonic analyses of the dimensionless 
modal response variables selected wherein QVJ(I), QV(K), and QT(J) are, 
respectively, the (I) flatwise, (K) edgewise and (J) torsional uncoupled 
mode responses. 

Sample Page 12 depicts the harmonic analyses of the total shears and 
moments exerted by one blade to the hub. In contrast to the steady hub loads 
listed in the AERODYNAMIC PERFORMANCE AND STRESSES output (Sample Page 9) 
which are calculated by integrating only the aerodynamic load distributions, 
the total hub loads which are herein harmonically analyzed are calculated by 
similarly integrating the combined aerodynamic and the dynamic load distribu- 
tions. The longitudinal, lateral and vertical hub shears comprising the first 
three quantities of this sample page all have the dimensions of lb and are 
defined in the x^-(aft), y^^- (starboard) , and zj^-(up and along axis of rota- 
tion) axis directions, respectively. The roll, pitch and yaw moments compris- 
ing the latter three quantities on this sample page have the dimensions of 
Ib-ft and are defined positive (using the right-hand rule) about the x^^-, 
y^^-, and zj^-axes, respectively. Note that the aerodynamic rolling moment 
whose output is depicted in Sample Page 9 is defined positive starboard side 
down and is opposite from the harmonically analyzed total rolling moment 
depicted in Sample Page 12. Sample Page 13 depicts the harmonic analysis of 
the flatwise stresses at the center of each of the spanwise segments, A 
similar output listing is provided for both edgewise and torsional stresses. 
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Impedance Matrices 


When the analysis Is optionally used to calculate Impedance matrices 
(as controlled by input items 930 through 939, and 946) these matrices are 
output as shown in Sample Pages 14 and 15. Sample Page 14 presents the output 
format for the blade passage harmonic characteristics of conventionally defined 
hub loads and their respective imepdances due to vibratory hub acceleration. 

The column denoted HUB LOADS presents the N/rev harmonic content of the hub 
loads as determined by summing contributions from all blades, where N is 
indicated in the first line of this sample page. Thus, the values in this 
column are equal to the Nth harmonic values shown in Sample Page 12 but 
multiplied by the number of blades. Note that, as shown in Sample Page 14, 
the results would be appropriate for the 4th harmonic hub loads for a four- 
bladed rotor. The columns denoted LDXC, LDYC, LDZC, LDXS, LDYS and LDZS 
present the derivatives of the respective HUB LOADS values with respect to 
N/rev vibratory displacements in the x, y and z directions, cosine and sine 
components each. The columns denoted RDXC, RDYC, ..., RDZS present similar 
derivatives of the hub loads with respect to vibratory rotations about the 
X, y and z axes, cosine and sine components each. Note that the units of 
the loads and perturbational quantities are defined in the second line of this 
sample page. The positive directions of the loads are also defined and are 
consistent with those of Sample Page 12. 

Sample Page 15 presents elements of an "augmented" impedance matrix 
representing the partial derivatives of the HUB LOADS components, but with 
respect to the six components of higher harmonic (pitch) control angle 
(HHC) appropriate to N/rev vibratory hub loads. Normally these components 
consist of the cosine and sine components each of the N-1, N and N+1 harmonics. 
Note that Sample Page 15 indicates harmonic components of HHC appropriate to 
4/rev vibrations. As indicated in the first line of this sample page the 
pitch angle perturbations have the units of degrees. 

In Sample Page 16 is depicted the typical additional page of output 
generated at the beginning of every case following the first case of a multiple 
case sum. The two columns depict, respectively, the location numbers and data 
values for the newly inputted data distinguishing the prsent case from the 
previous one. This feature is intended solely as an ease of usage output to 
assist in data management. 
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SUMMARY OF ROTOR HUB A/REV LOADS AND HUB IMPEDANCE MATRIX CALCULATIONS FOR A/REV VIBRATORY HUB DISPLACEMENTS 
UNITS... U(SHEAR) - LBS, U(MOMENT) - FT-LBS, U(LINEAR DISPLACEMENT) « FT, U(ROTATIONAL DISPLACEMENT) - RAD 
CONVENTIONAL FOURIER SERIES REPRESENTATIONS 







HUB LOADS 

LDXC 

LDYC 

LDZC 

RDXC 

RDYC 

RDZC 

LONG. 

SHR 

(+» 

AFT) 
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X.XXXE-YY 

X.XXXE-YY 
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X.XXXE-YY 
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X.XXXE-YY 

LAT. 

SHR 

(+, 

STRBD) 

SIN 

X.XXXE-YY 

X.XXXE-YY 
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HUB IMPEDANCE MATRIX FOR HHC PITCH ANGLES (UNITS - DEG) 
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X.XXXE-YY 
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APPENDIX I 


UPDATES TO DOCUMENTATION OF BASIC G400 ANALYSIS 

Since the publication of Reference 1 the original G400 analysis has 
undergone continual development. In this process many of the original forimi- 
lations and derivations have been reassessed, recheckcd , and, where appropr J/ite , 
reformulated. This appendix presents some of these reformulations and are, 
in fact, an errata update sheet to Reference 1. 

1. The first such update pertains to the definitions for the deflection 
correction terms resulting from structural twist. In particular, the first 
order terms given by Eqs. (7) and (8) of Reference 1 should read, respectively, 
as: 



(I.l) 


( 1 . 2 ) 


2. The second update pertains to formulation of the elements of the 
inertia coupling matrix, [A], (Eq. (45) of Reference 1). This matrix should 

3nd 


be completely symmetric; the offending terms, and A^, should equal 


their respective transposed elements: 


*“i»i 




(1.3) 






AW|^ dr 


(1.4) 


3. The third and fourth updates pertain to the implementation of the 
redundant analysis of flexbeam-torque tube assemblies. The first of these two 
relates to the torsional stiffness characteristics of the flexbeams. A more 
correct statement of the boundary conditions for the flexbeam torsion equation 
(Eq. (115) of Reference 1) includes relating the flexbeam outboard twist rate, 
0'(r.), to the twist rate of the outer blade: 


e'(r^)-©' 


(1.5) 
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where the twist rate of the outer blade comes from built-in and response 
contributions: 

NTM 

©j=©B(rj)+Z (i.6) 

js| J 

The inclusion of this generally nonzero value of twist rate has the effect 
of modifying the resulting expression for the elastic torsion deflection of the 
flexbeam at the juncture (Eq. (129a)); 


^®J ' (61,623-6216,3) ^ [®I3^^2 * i+| ®22 ®J^ 

-e23<fl.j+, +e,2©')] 


( I .7) 


and the symbolic expression for this deflection (Eq. (130b)) then becomes: 


+ + (1-8) 

+ 00'. ©j 

4. The fourth update pertains to the calculation of internal blade shears 
and moments in the blade adjacent to the juncture. As outlined in Reference 1 
a force integration approach was implemented which necessitated extrapolations 
to estimate the modal accelerations. This calculation was found to be best 
performed using a mode deflection method Instead. Thus, Eqs. (132b) and 132f). 
of Reference 1 are, respectively, replaced by the following expressions: 

Sy5B="S®J^, DSI, q„j+C0jI^ DS3kQvj^ (1. 9a) 


NEM 

Sz5^-*c€^jI DSI.q^, +s0jZ DS3nqv. (1. 9b) 

6 Jr| ' kM 


NTM 
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where: 



NFM 

= C©j I DS2|qwi-S® Z 
B ''isi ' ‘'ksi 


DS4r Qv^ 


NFM NEM 

Mz5b= 


DS\.^--ul f my dr 

DS2, --J m[{7--i^)wJ,.y„.-T{y^.(r)-y„.(rj))]d7 
•■j 



dr 


DS4jj=J’ m ^ " Xv,^t dr 

Tj 



(i.9d) 


(I.9e) 


(I.lOn) 


(I. 10b) 


(1.10c) 


(I.lOd) 


(I.lOe) 
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